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first  two  features)  and  moderate  to  high  thrust  (limited  primarily  by  the 
average  laser  power  available ) •  The  present  research  addresses  questions 
related  to  thruster  performance  and  optical  design.  j 

In  the  thruster  scheme  under  consideration,  parabolic  nozsle  walls  focus 
the  incoming  laser  beam  to  yield  breakdown  in  a  propellant  at  the  focal 
point  of  the  parabola.  The  resulting  high  pressure  plasma  is  character is* 
tic  of  a  detonation  wave  initiation  by  high  power  laser-induced  breakdown. 
With  a  short  laser  pulse,  the  detonation  wave  quickly  becomes  a  blast  wave 
which  propagates  to  the  nozzle  exit  plane  converting  the  high  pressure  of 
the  gas  behind  it  to  a  force  on  the  nozzle  wall.  Propellant  is  fed  to  the 
focal  region  from  a  plenum  chamber.  The  laser-induced  blast  wave  stops 
the  propellant  flow  through  the  throat  until  the  pressure  at  the  throat 
decays  to  the  sonic  pressure;  then  the  propellant  flow  restarts.  The 
process  is  repeated  with  each  successive  laser  pulse. 

j 

A  major  part  of  the  present  work  involved  carrying  out  thruster  perfor¬ 
mance  experiments  using  a  small  scale  conical  nozzle,  external  focusing 
optics,  and  pulsed  CO2  (X-10.6  ym)  and  pulsed  XeF  (X-0.35  urn)  lasers  as 
the  propellant  heating  sources.  All  experiments  were  performed  in  a 
vacuum  environment.  With  the  XeF  laser  (Ep*2.7  J,  Tp-0.6  ms)  and  argon 
propellant,  a  maximum  specific  impulse  of  500  seconds  was  achieved  with  a 
conversion  efficiency  of  laser  energy  to  propellant  kinetic  energy  upwards 
of  50%.  With  the  same  propellant  and  the  CO2  TEA  laser  (Ep“6.5  J,  Tp*3  y 
ys),  the  maximum  specific  impulse  was  700  seconds  with  a  laser  to  kinetic 
energy  conversion  efficiency  "  20%.  Finally,  using  hydrogen  propellant 
and  the  CO2  laser  as  the  heating  source,  the  experimental  data  indicate 
specific  impulse  up  to  3000  seconds  with  corresponding  energy  conversion 
efficiencies  upwards  of  *  70%. 

.4 

•  f 

In  -support  of  the  thruster  performance  experiments,  a  detailed  computer 
mo^fl  that  was  developed  previously  was  modified  and  used  to  simulate  the 
-no^le  fluid  dynamics  in  the  conical  thruster.  The  model  treats  non- 
st^dy  nozzle  flow  of  a  real  gas  with  laser  heat  addition.  Major  enphasis 
in  the  present  work  was  directed  toward  simulating  the  fluid  dynamics  ap¬ 
propriate  for  repetitively-pulsed  operation,  i.e.,  the  situation  in  which 
the  nozzle  has  only  partly  refilled  with  fresh  gas  between  successive 
laser  pulses  and  the  remaining  downstream  section  of  the  nozzle  contains 
the  low  pressure  tail  of  the  previous  blast  wave.  For  a  given  initial 
propellant  mass  distribution  and  laser  input,  the  computer  program  calcu¬ 
lates  such  quantities  as  the  exhaust  gas  momentum,  kinetic  energy,  and 
specif ic  impulse . 

Finally,  an  assessment  has  been  made  of  optical  design  criteria  for  a 
self-focusing  paraboloidal  nozzle*  The  primary  conclusions  reached  are 1 
( 1 )  the  mission  requirements  for  off-axis  performance  allow  the  surface  to 
be  of  only  very  modest  optical  quality,  (2)  for  off-axis  angles  i  10°, 
only  very  modest  magnifications  (S  25)  can  be  anticipated,  (3)  optimum  op¬ 
tical  and  fluid  dynamic  performance  appear  to  be  mutually  Inclusive,  con¬ 
sequently,  a  compromise  is  required,  and  (4)  the  impact  of  off-axis  focus¬ 
ing  and  breakdown  on  thruster  performance  should  be  evaluated. 
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In  this  report  we  present  results  of  continuing  analytical  and  experi¬ 
mental  investigations  carried  out  to  evaluate  the  concept  of  pulsed  laser  pro¬ 
pulsion.  This  advanced  propulsion  scheme,  which  has  been  the  subject  of  sev¬ 
eral  previous  studies,  involves  supplying  propellant  energy  by  beaming  short, 
repetitive  laser  pulses  to  a  thruster  from  a  remote  laser  power  station.  The 
concept  offers  the  advantages  of  a  remote  power  source,  high  specific  isqmlse, 
high  payload  to  total  mass  ratio  (a  consequence  of  the  first  two  features)  and 
moderate  to  high  thrust  (limited  primarily  by  the  average  laser  power  avail¬ 
able).  The  present  research  addresses  questions  related  to  thruster  perfor¬ 
mance  and  optical  design. 

In  the  thruster  scheme  under  consideration,  parabolic  nozzle  walls 
focus  the  incoming  laser  beam  to  yield  breakdown  in  a  propellant  at  the  focal 
point  of  the  parabola.  The  resulting  high  pressure  plasma  is  characteristic 
of  a  detonation  wave  initiation  by  high  power  laser-induced  breakdown.  With  a 
short  laser  pulse,  the  detonation  wave  quickly  becomes  a  blast  wave  which 
propagates  to  the  nozzle  exit  plane  converting  the  high  pressure  of  the  gas 
behind  it  to  a  force  on  the  nozzle  wall.  Propellant  is  fed  to  the  focal 
region  from  a  plenum  chamber.  The  laser-induced  blast  wave  stops  the  propel¬ 
lant  flow  through  the  throat  until  the  pressure  at  the  throat  decays  to  the 
sonic  pressure;  then  the  propellant  flow  restarts.  The  process  is  repeated 
with  each  successive  laser  pulse. 

A  major  part  of  the  present  work  involved  carrying  out  thruster  per¬ 
formance  experiments  using  a  small  scale  conical  nozzle,  external  focusing  op¬ 
tics,  and  pulsed  CO2  (X-10.6  urn)  and  pulsed  XeF  (1-0.35  um)  lasers  as  the  pro¬ 
pellant  heating  sources.  All  experiments  were  performed  in  a  vacuum  environ¬ 
ment.  With  the  XeF  laser  (Ep«2.7  J,  Tp-0.6  us)  and  argon  propellant,  a  maxi¬ 
mum  specific  impulse  of  500  seconds  was  achieved  with  a  conversion  efficiency 
of  laser  energy  to  propellant  kinetic  energy  upwards  of  50%.  With  the  same 
propellant  and  the  CO2  TEA  laser  (Ep"6.5  J,  Tp-3  us),  the  maximum  specific 
impulse  was  700  seconds  with  a  laser  to  kinetic  energy  conversion  efficiency 
*  20%.  Finally,  using  hydrogen  propellant  and  the  CX>2  laser  as  the  heating 
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source,  the  experimental  data  indicate  specific  impulse  up  to  3000  seconds 
with  corresponding  energy  conversion  efficiencies  upwards  of  ■  70%. 

In  support  of  the  thruster  performance  experiments,  a  detailed  comput¬ 
er  model  that  was  developed  previously  was  modified  and  used  to  simulate  the 
nozzle  fluid  dynamics  in  the  conical  thruster.  The  model  treats  non-steady 
nozzle  flow  of  a  real  gas  with  laser  heat  addition.  Major  emphasis  in  the 
present  work  was  directed  toward  simulating  the  fluid  dynamics  appropriate  for 
repetitively-pulsed  operation,  i.e.,  the  situation  in  which  the  nozzle  has 
only  partly  refilled  with  fresh  gas  between  successive  laser  pulses  and  the 
remaining  downstream  section  of  the  nozzle  contains  the  low  pressure  tail  of 
the  previous  blast  wave.  For  a  given  initial  propellant  mass  distribution  and 
laser  input,  the  computer  program  calculates  such  quantities  as  the  exhaust 
gas  momentum,  kinetic  energy,  and  specific  impulse. 

Finally,  an  assessment  has  been  made  of  optical  design  criteria  for  a 
self-focusing  paraboloidal  nozzle.  The  primary  conclusions  reached  are: 

(1)  the  mission  requirements  for  off-axis  performance  allow  the  surface  to  be 
of  only  very  modest  optical  quality,  (2)  for  off-axis  angles  i  10°,  only  very 
modest  magnifications  (S  25}  can  be  anticipated,  (3)  optimum  optical  and  fluid 
dynamic  performance  appear  to  be  mutually  inclusive,  consequently,  a  compro¬ 
mise  is  required,  and  (4)  the  impact  of  off-axis  focusing  and  breakdown  on 
thruster  performance  should  be  evaluated. 
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Schematic  of  repetitively  pulsed  laser-powered  thruster. 

2.1  Density  at  2  cm  in  real  argon,  for  stagnation  conditions  11 

of  300  K  and  10  atm,  with  1  J  deposited  between  the  throat 
and  0.1  cm.  Density  reduced  by  factors  of  1,  100  and  1000 
at  2  cm  before  starting  the  computation. 

2.2a  Mass  averaged  velocity  distributions  at  various  times.  No  13 

density  reduction.  Real  argon  for  stagnation  conditions  of 
300  K  and  10  atm,  with  1  J  deposited  between  the  throat  and 
0. 1  cm. 

2.2b  Mass  averaged  velocity  distributions  at  various  times.  14 

Density  reduced  by  a  factor  of  1000  downstream  of  2  cm. 

Real  argon  for  stagnation  conditions  of  300  K  and  10  atm, 
with  1  J  deposited  between  the  throat  and  0.1  cm. 

2.3a  Mass  fluence  distribution  at  various  times.  No  density  15 

reduction.  Real  argon  for  stagnation  conditions  of  300  K 
and  10  atm,  with  1  J  deposited  between  the  throat  and  0.1 
cm. 

2.3b  Mass  fluence  distributions  at  various  times.  Density  16 

reduced  by  1000  at  downstream  of  2  cm.  Real  argon  for 
stagnation  conditions  of  300  K  and  10  atm,  with  1  J 
deposited  between  the  throat  and  0.1  cm. 

2.4a  Momentum  fluence  distributions  at  various  times.  No  density  18 

reduction.  Real  argon  for  stagnation  conditions  of  300  K 
and  10  atm,  with  1  J  deposited  between  the  throat  and  0.1 
cm. 

2.4b  Momentum  fluence  distributions  at  various  times.  Density  19 

reduced  by  a  factor  of  1000  downstream  of  2  cm.  Real  argon 
for  stagnation  conditions  of  300  K  and  10  atm,  with  1  J 
deposited  between  the  throat  and  0.1  cm. 

2.5a  Gas  velocity  distributions  at  various  times.  No  density  20 

reduction.  Real  argon  for  stagnation  conditions  of  300  K 
and  10  atm,  with  1  J  deposited  between  the  throat  and  0.1 
cm* 

2.5b  Gas  velocity  distributions  at  various  times.  Density  21 

reduced  by  a  factor  of  1000  downstream  of  2  cm.  Real  argon 
for  stagnation  conditions  of  300  K  and  10  atm,  with  1  J 
deposited  between  the  throat  and  0. 1  cm. 
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2.6a  Pressure  history  at  pressure  transducer  locations.  No  23 

density  reduction.  Real  argon  for  stagnation  conditions  of 
300  K  and  10  atm,  with  1  J  deposited  between  the  throat  and 
0. 1  cm. 

2.6b  Pressure  history  at  pressure  transducer  locations.  Density  24 

reduced  by  a  factor  of  1000  downstream  of  2  cm.  Real  argon 
for  stagnation  conditions  of  300  K  and  10  atm.  with  1  J 
deposited  between  the  throat  and  0.1  cm. 

2.7  Shock  trajectory  in  quiescent  air.  expressed  in  similarity  26 

variables.  Eq.  (2.2).  Perfect  air,  y-1.4,  with  ambient 
conditions  of  1  atm  and  300  K.  Hot  gas  initial  state  4000 

K.  ambiert  density,  except  for  squares,  where  initial  state 
is  ambient  density,  121  atm  pressure. 

2.8  Peak  shock  pressure  as  a  function  of  distance,  expressed  in  28 

similarity  variables,  Eqs.  (2.2)  and  (2.5).  Perfect  air, 

Y*1.4,  with  ambient  conditions  of  1  atm  and  300  K.  Hot  gas 
initial  state  4000  K,  ambient  density,  except  for  squares, 
where  initial  state  is  ambient  density,  1 21  atm  pressure. 

3.1  Schematic  diagram  of  laser  propulsion  experiment:  32 

C  -  Luminics  SOD  calorimeter;  QW  -  quartz  window; 

M  -  aluminum  coated  concave  surface  mirror  (f-160  cm); 

OAP  -  off-axis  parabola  concave  mirror  (f*15  cm); 

RC  -  rocket  cone;  P  -  plenum,  G  -  gold  surface  mirror; 

KC1  -  KC1  window;  PZ  -  Piezo  pressure  transducer; 

SV  -  Solenoid  valve. 

3.2  Schematic  of  nozzle  configuration  for  laser  propulsion  35 

experiment. 

3.3  Representative  Oscilloscope  records  of  the  XeF  laser  energy  40 

measurements  (-a),  the  pressure  history  in  the  plenum  (-b), 

and  the  laser-induced  shock  wave  measurements  (-c).  (a)  50 

ms/div,  50mV/div:  (b)  10  ms/div  100  mV/div;  discontinuity  on 
the  plateau  region  indicates  the  time  when  the  CO2  and  the 
XeF  lasers  were  fired:  (c)  10  ws/div,  upper  500  mV,  lower  100 
mV/div,  respectively;  +  indicates  the  CO2  laser  firing;  *  and 
•  indicate  the  arrival  of  the  CO2  laser-induced  blast  wave  at 
the  first  and  second  pressure  transducer  station,  respective¬ 
ly;  indicates  the  XeF  laser  firing;  A  and  £  indicate  the 
arrival  of  the  XeF  laser-induced  shock  wave  at  the  first  and 
second  pressure  transducer  station,  respectively. 
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3.4  Blast  wave  measurement  In  1  atm  air  using  the  CO2  TEA  laser  44 

(a)  and  the  XeF  laser  (b).  (a)  Pressure  profiles  at  the 

first  and  the  second  pressure  transducer  stations,  respect¬ 
ively.  Beginning  of  oscilloscope  trace  corresponds  to 
firing  of  the  laser:  20  ps/div:  upper  trace  -  2  V/div;  lower 
trace  -  500  mV/div.  (b)  The  XeF  laser-induced  blast  wave 
measurement:  20  ps/div j  upper  trace  -  2  V/div;  lower  trace  - 
500  mV/div. 

3.5  Blast  wave  energy  measurement  in  argon  using  the  XeF  laser.  46 

(a)  XeF  laser  output  energy  measurement:  10  ms/div;  50 

mV/div.  (b)  Transit  time  measurement  of  the  XeF  laser-induced 
blast  wave;  5  ps/div;  upper  trace  -  100  mV/div;  lower 
trace  -  50  mV/div,  *  and  •  indicate  the  arrival  of  blast  wave 
at  the  first  and  second  stations,  respectively. 

3.6  Conversion  efficiency  of  the  CO2  laser  (10.6  pm)  energy  to  47 

blast  wave  energy  versus  argon  mass  flow  rate  through  the 

nozzle.  Error  bar  indicates  the  maximum  error  incurred  in 
experiments,  evaluated  using  equation  (3.4.8). 

3.7  Conversion  efficiency  of  the  XeF  laser  (0.35  pm)  energy  to  48 

blast  wave  energy  versus  argon  flow  rate  through  the  nozzle. 

Error  bar  indicates  the  maximum  error  incurred  in  experiments, 
evaluated  using  Eq.  (3.4.8). 

3.8  Gas  kinetic  energy  conversion  efficiency  versus  mass  with  50 

argon  mass  flow  rate  of  2.22  g/s.  Dashed  lines  correspond  to 

ISp(s)  “  500  ( - ),  400  ( —  — ),  350  ( —  •  — )  with  XeF  energy 

of  2.6  J.  Error  bar  indicates  the  maximum  error  incurred  in 
experiment,  which  was  evaluated  using  Eq.  (3.4.9). 

3.9  Gas  kinetic  energy  conversion  efficiency  versus  mass  with  51 

argon  mass  flow  rate  of  1.89  g/s.  Dashed  lines  correspond  to 

Isp(s)  -  500  ( - ),  400  ( - ),  350  ( —  •  — )  with  the  XeF 

energy  of  2.6  J.  Error  bars  indicate  the  maximum  error 

incurred  in  experiment,  which  was  evaluated  using  Eq.  (3.4.9). 

3.10  Gas  kinetic  energy  conversion  efficiency  versus  mass  with  52 

argon  mass  flow  rate  of  1.55  g/s.  Dashed  lines  correspond  to 

I8p(s)  "  500  ( - ),  400  ( - ),  350  (—  •  — )  with  the  XeF 

energy  of  2.6  J.  Error  bar  indicates  the  maximum  error 

incurred  in  experiments,  which  was  evaluated  using  Eq. 

(3.4.9). 
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3.11  Gas  kinetic  energy  conversion  efficiency  versus  mass  with 
argon  mass  flow  rate  of  1.32  g/s.  Dashed  lines  correspond 

to  Isp(s)  **  500  ( - ),  400  ( - ),  350  (—  •  — •)  with  Xef 

energy  of  2.6  J.  Error  bar  indicates  the  maximum  error 
incurred  in  experiments ,  which  was  evaluated  using  Eq. 
(3.4.9). 

3.12  Gas  kinetic  energy  conversion  efficiency  versus  mass  with 
argon  mass  flow  rate  of  0.96  g/s.  Dashed  lines  correspond 

to  Isp(s)  ■  500  ( - ),  400  ( —  — ),  350  ( —  •  — )  with  the 

XeF  energy  of  2.6  J.  Error  bar  indicates  the  maximum  error 
incurred  in  experiments ,  which  was  evaluated  using  Eq. 
(3.4.9). 

3.13  Summary  of  efficiencies  achieved  for  conversion  of  XeF  laser 
energy  into  kinetic  energy  of  argon  propellant.  Results  are 
plotted  as  a  function  of  propellant  mass  flow  rate.  Error 
bars  correspond  to  1  a  of  standard  deviation. 

3.14  Specific  impulse  data  obtained  for  argon  propellant  heated 
by  two  C02  TEA  laser  pulses.  Specific  impulse  is  determined 
from  transit  time  measurements  of  pressure  wave  induced  by 
2nd  laser  pulse.  Dashed  lines  plotted  with  data  indicate 
loci  of  constant  propellant  kinetic  energy  (1/2  Am(g  Isp) 2 ) . 
Am  is  the  mass  expelled  through  throat  between  first  and 
second  laser  pulses. 

3.15  Specific  impulse  data  obtained  for  hydrogen  propellant 
heated  by  two  C02  TEA  laser  pulses.  Specific  impulse  is 
determined  from  transit  time  measurements  of  pressure  wave 
induced  by  second  laser  pulse.  Dashed  lines  are  loci  of 
constant  propellant  kinetic  energy  (1/2  Am(g  Isp)2).  Am  is 
the  mass  expelled  through  throat  between  first  and  second 
laser  pulses. 
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INTRODUCTION  AND  EXECUTIVE  SUMMARY 


Repetitively-pulsed  (RP)  laser  propulsion  Is  a  pulsed  jet  propulsion 
concept  in  which  the  propellant  energy  is  supplied  by  the  absorption  of  short, 
repetitive  laser  pulses  beamed  to  the  thruster  from  a  remote  laser  power  sta¬ 
tion.  This  advanced  propulsion  concept  offers  the  advantages  of  a  remote 
power  source,  high  specific  impulse,  high  payload  to  total  mass  ratio  (a  con¬ 
sequence  of  the  first  two  features),  and  moderate  to  high  thrust  (limited  pri¬ 
marily  by  the  average  laser  power  available ) . 

In  the  RP  thruster  concept  shown  in  Fig.  1.1,  parabolic  nozzle  walls 
focus  the  incoming  beam  to  yield  breakdown  in  a  propellant  at  the  focal  point 
of  the  parabola.  The  resulting  high  pressure  plasma  is  characteristic  of  det¬ 
onation  wave  intiation  by  high  power  laser-induced  breakdown.  With  a  short 
laser  pulse,  the  detonation  wave  quickly  becomes  a  blast  wave,  which  propa¬ 
gates  to  the  nozzle  exit  plane  converting  the  high  pressure  of  the  propellant 
gas  behind  it  into  a  force  on  the  nozzle  wall.  Propellant  is  fed  to  the  focal 
region  from  a  plenum  chamber .  The  laser-induced  blast  wave  stops  the  propel¬ 
lant  flow  through  the  throat  until  the  pressure  at  the  throat  decays  to  the 
sonic  pressure;  then,  the  propellant  flow  restarts.  The  process  is  repeated 
with  each  successive  laser  pulse. 

Physical  Sciences  Inc.  (PSI)  has  been  developing  the  technology  of 
pulsed  laser  propulsion  under  DARPA  support.1,1”1,4  Previous  experiments  and 
supporting  theoretical  studies  have  shown  that  repetitively  pulsed  10.6  urn 
laser  radiation  can  be  utilized  to  heat  simple  propellants  to  high  stagnation 
temperatures;  the  result  is  rocket  propulsion  at  high  specific  impulse  (ISp). 
Using  a  simple  thruster  concept  (similar  to  the  one  shown  in  Fig.  1.1)  with  a 
variety  of  propellants,  we  obtained  specific  impulses  of  up  to  1000  s  with  an 
energy  conversion  efficiency  of  40%.  PSI  has  since  been  engaged  in  utilizing 
the  knowledge  and  experience  obtained  from  this  program  to  advance  laser  pro¬ 
pulsion  technology  for  DARPA  space  applications.  Since  laser  propulsion  ac¬ 
tivities  are  driven  by  advances  in  laser  technology  for  other  applications  in 
space,  it  has  been  proposed  that  pulsed  laser  propulsion  systems  utilize  vis- 
ible/UV  pulsed  lasers.  Therefore,  the  recent  program  objectives  have  been  to 
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1)  establish  missions  and  associated  laser  requirements  for  a  ground  based 
laser  propulsion  system,  2)  establish  a  thruster  design  capability  for  effi¬ 
cient  performance  as  a  function  of  desired  specific  impulse  and  laser  wave¬ 
length,  3)  develop  scaling  laws  for  laser  propulsion  with  pulsed  visible 
lasers,  and  4)  identify  the  technology  issues  for  laser  propulsion  with  vis¬ 
ible  lasers. 

The  technical  approach  for  the  previous  program  phase14  involved  per¬ 
forming  theoretical  modeling  and  systems  analysis  along  with  supporting  exper¬ 
imental  studies  using  a  commercially  built  pulsed  XeF  laser.  Although  it  has 
not  been  established  that  XeF  (0.35  urn)  is  indeed  the  desirable  wavelength, 
the  advanced  nature  of  excimer  laser  technology  has  led  to  the  availability  of 
commercially  built  lasers  with  sufficient  energy  and  the  proper  pulse  time  to 
perform  meaningful  scaling  experiments.  The  analytical  aspects  of  this  pro¬ 
gram  were  divided  into  studies  of  propulsion  physics  wavelength  scaling  and 
mission  analysis  to  identify  DOD  applications  of  laser  propulsion  systems. 

The  propulsion  physics  studies  were  carried  out  at  PSI,  and  the  mission  analy¬ 
sis  was  performed  both  at  PSI  and  the  Lockheed  Missiles  and  Space  Company 
under  subcontract  to  PSI.  The  experimental  program  was  performed  at  PSI  using 
in-house  pulsed  lasers. 

As  a  result  of  the  mission  analysis,  several  attractive  missions  were 
identified  with  perhaps  the  most  attractive  near  term  one  being  the  orbital 
repositioning  of  satellites.  In  addition,  the  absorption  physics  scaling  ex¬ 
periments  served  to  establish  the  wavelength  scaling  of  breakdown  and  absorp¬ 
tion  in  several  gases  of  interest.  For  argon,  the  experimental  measurements 
were  also  supported  by  corresponding  theoretical  model  calculations  of  the 
plasma  formation  process.  Finally,  preliminary  single  pulse  0.35  pm  thruster 
performance  data  were  acquired  and  compared  with  numerical  code  simulations  of 
the  experiments  that  model  the  laser  absorption  physics,  equilibrium  gas  chem¬ 
istry  (including  multiple  ionization)  and  gas  dynamics  expansion.  The  latter 
data  indicated  that  laser  to  blast  wave  energy  conversion  efficiencies  in  ex¬ 
cess  of  50%  are  possible  at  this  wavelength. 

For  the  present  program,  work  was  carried  out  in  three  areas:  1)  ex¬ 
perimental  thruster  performance  experiments  to  extend  the  previous  single- 
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pulse  0.35  im  data  to  multiple  pulse  operation,  2)  modification  of  the  numer¬ 
ical  flow  code  to  simulate  the  fluid  dynamics  appropriate  to  repetitively- 
pulsed  operation,  i.e.,  a  nozzle  that  has  only  partly  refilled  with  fresh  gas 
prior  to  each  laser  pulse,  and  3)  an  analytical  assessment  of  the  optical  per¬ 
formance  of  a  self-focusing  paraboloidal  nozzle. 

The  thruster  performance  experiments  were  carried  out  using  a  small 
scale  conical  nozzle,  external  focusing  optics,  and  pulsed  CO2  (X»10.6  urn)  and 
pulsed  XeF  (X*0.35  um)  lasers  as  the  propellant  heating  sources.  All  experi¬ 
ments  were  performed  in  a  vacuum  environment.  With  the  XeF  laser  (Ep*2.7  J, 
Tp«0.6  us)  and  argon  propellant,  a  maximum  specific  impulse  of  500  seconds  was 
achieved  with  a  conversion  efficiency  of  laser  energy  to  propellant  kinetic 
energy  upwards  of  50%.  With  the  same  propellant  and  the  C02  TEA  laser 
(Ep"6. 5  J,  Tp»3  us),  the  maximum  specific  impulse  was  700  seconds  with  a  laser 
to  kinetic  energy  conversion  efficiency  *  20%.  Finally,  using  hydrogen  pro¬ 
pellant  and  the  C02  laser  as  the  heating  source,  the  experimental  data  indi¬ 
cate  specific  impulse  up  to  3000  seconds  with  corresponding  energy  conversion 
efficiencies  upwards  of  "  70%.  The  experimental  results  are  presented  in  Sec¬ 
tion  3. 


In  support  of  the  thruster  performance  experiments,  the  detailed  com¬ 
puter  model  that  was  developed  previously  was  modified  and  used  to  simulate 
the  nozzle  fluid  dynamics  in  the  conical  thruster.  The  model  treats  non¬ 
steady  nozzle  flow  of  a  real  gas  with  laser  heat  addition.  Major  emphasis  in 
the  present  work  was  directed  toward  simulating  the  fluid  dynamics  appropriate 
for  repetitively-pulsed  operation,  i.e.,  the  situation  in  which  the  nozzle  has 
only  partly  refilled  with  fresh  gas  between  successive  laser  pulses  and  the 
remaining  downstream  section  of  the  nozzle  contains  the  low  pressure  tail  of 
the  previous  blast  wave.  For  a  given  initial  propellant  mass  distribution  and 
laser  input,  the  computer  program  calculates  such  quantities  as  the  exhaust 
gas  momentum,  kinetic  energy,  and  specific  impulse.  The  model  calculations 
are  described  in  Section  2. 

Finally,  an  assessment  has  been  made  of  optical  design  criteria  for  a 
self -focusing  paraboloidal  nozzle.  The  primary  conclusions  reached  are: 

1)  the  mission  requirements  for  off-axis  performance  allow  the  surface  to  be 
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of  only  very  modest  optical  quality,  2)  for  off-axis  angles  l  10°,  only  very 
modest  magnifications  <$  25)  can  be  anticipated,  3)  optimum  optical  and  fluid 
dynamic  performance  appear  to  be  mutually  exclusive,  consequently  a  compromise 
is  required,  and  4)  the  impact  of  off-axis  focusing  and  breakdown  on  thruster 
performance  should  be  evaluated.  The  details  of  the  optical  assessment  task 
are  presented  in  Section  4. 
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2.  MODELING  OF  A  PULSED  LASER-HEATED  THRUSTER  -  SINGLE  PULSE 
AND  MULTIPLE  PULSE  OPERATION 

2. 1  Introduction 

The  main  purpose  of  the  modeling  effort  was  to  calculate  the  flow  in  a 
nozzle  which  was  only  partly  filled  with  fresh  gas.  This  simulates  the  case  in 
which  the  time  between  laser  pulses  is  short  enough  so  that  the  nozzle  does  not 
refill  completely  after  the  flow  from  the  stagnation  chamber  restarts.  The  re¬ 
gion  of  the  nozzle  near  the  throat  contains  the  fresh  gas  which  has  entered 
through  the  throat.  However,  the  rest  of  the  nozzle  contains  heated  gas  left 
over  from  the  previous  pulse,  which  is  at  much  lower  density. 

An  estimate  of  the  time  necessary  to  fill  the  nozzle  can  be  obtained 
from  the  ratio  of  the  nozzle  length  to  the  limiting  velocity  of  the  incoming 
gas.  For  300  K  argon  the  limiting  velocity  is  5.59  x  io4  cm/s,  so  for  a  10  cm 
nozzle  length  it  takes  180  us  to  fill  the  nozzle.  The  time  for  the  laser¬ 
generated  blast  wave  to  traverse  10  cm  might  be  another  10  us.  So  if  the  inter¬ 
pulse  time  is  less  than  200  us  (more  than  5000  pulses  per  second)  the  tail  of 
the  previous  pulse  will  still  be  in  a  10  cm  nozzle  when  the  next  pulse  occurs. 
Such  conditions  were  explored  in  the  experiments  ( see  Section  3 ) . 

The  purpose  of  having  a  pulse  repetition  rate  higher  than  the  rate 
which  permits  complete  refilling  is  to  increase  the  average  specific  impulse. 

It  is  clear  that  the  fluid  mechanics  will  be  different  if  part  of  the  nozzle  is 
filled  with  the  lower  density  gas  representative  of  the  tail  of  the  previous 
pulse,  instead  of  the  fresh  gas  which  has  come  in  through  the  throat. 

An  early  study  of  the  fluid  mechanics  of  pulsed  laser  propulsion  by 
Simons  and  Pirri  (Ref.  2.1)  dealt  with  the  multi-pulse  problem.  However,  the 
valuable  analytical  modelling  in  that  paper  was  restricted  to  a  perfect  gas,  and 
utilized  blast  wave  theory,  in  which  the  shock  wave  is  always  hypersonic 
relative  to  the  gas  in  the  nozzle.  In  the  work  described  in  Ref.  2.2  these  two 
restrictions  were  removed.  A  computer  model  of  the  non-steady  nozzle  flow  of  a 
real  gas  with  laser  heat  addition  was  constructed,  and  used  to  perform  single 
pulse  calculations. 
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During  the  course  of  the  present  work,  this  computer  program,  described 
in  Ref.  2.2,  was  modified  to  simulate  the  flow  in  the  partially-filled  nozzle 
and  some  runs  were  made  to  compare  the  filled  and  partially  filled  cases,  ttiis 
part  of  the  work  is  described  in  the  next  section. 

The  program  was  also  used  to  aid  the  experimental  diagnostics  by  pro¬ 
viding  shock  wave  position  vs.  time  (trajectories)  as  well  as  shock  peak  pres¬ 
sure  vs.  shock  position,  for  some  early  experiments  in  quiescent  air.  The  cor¬ 
relations  resulting  from  these  calculations  are  presented  in  the  second  section 
below. 


2.2  Second  Pulse  Simulation 

Some  additional  output  of  the  LSDNS  computer  program  was  developed  in 
order  to  aid  the  interpretation  of  the  results  of  the  calculations.  To  motivate 
the  additional  output,  the  quasi -one -dimensional  conservation  equations  which 
are  solved  are  restated  here: 


Mass : 


Momentum : 


Energy: 


+  3PUA  ,  0 
3t  3x 


3puA  ,  3  , „ ,  ,  2»,  dA 

IT  +  3^  [Mp+pu  )]  ~ 


[pA(e+u2/2)]  +  [puA(h+u2/2)] 


0 


The  equations  are  written  in  conservation  form.  The  quantities  in  the  time  de¬ 
rivative  terms  are  usually  called  the  conserved  quantities,  while  those  in  the 
space  derivative  terms  are  the  fluxes  of  mass,  momentum  and  total  enthalpy. 


The  computer  program  was  modified  to  print  and  plot  the  x-distributions 
of  pA,  puA,  pA(e+u2/2),  and  the  x-distributions  of  the  time  integrals  of  puA, 
pAu2,  and  puA(h+u2/2) .  The  latter  three  flux  integrals  might  be  called  the 
mass,  momentum  and  energy  fluences. 


The  average  exhaust  velocity  for  a  pulsed  propulsion  device  can  be  de¬ 
fined  in  terms  of  a  mass-averaged  exhaust  velocity  and  the  mass  fluence  as 


u 

avg 
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where  the  integrals  are  evaluated  at  the  exhaust  plane.  The  denominator  is  the  ex¬ 
pelled  mass,  the  numerator  the  exhaust  momentum. 

To  simulate  the  presence  of  the  tail  of  the  previous  pulse  in  the 
nozzle,  a  reduced  density  region  is  required.  The  computer  program  initially 
fills  the  nozzle  with  steady  supersonic  flow  from  given  stagnation  conditions, 
assuming  a  perfect  gas.  A  simple  way  to  introduce  a  reduced  density  region  into 
this  flow  is  to  merely  divide  the  density  by  some  given  factor  in  the  region 
downstream  of  the  desired  nozzle  station.  However,  it  is  not  desirable  to  in¬ 
troduce  a  pressure  discontinuity,  since  the  calculation  scheme  will  treat  that 
as  an  interface  to  be  resolved  by  another  shock  wave,  which  is  not  the  intent. 

To  avoid  a  pressure  discontinuity  wherever  the  density  is  reduced,  the  internal 
energy  can  be  increased  by  the  same  factor  (equivalent  for  a  perfect  gas  to  in¬ 
creasing  the  temperature).  This  will  keep  the  variation  of  pressure  continu¬ 
ous. 


This  idea  was  implemented  in  the  computer  program  by  introducing  as  new 
input  quantities  a  station  at  which  the  density  and  internal  energy  change,  and 
the  factor  by  which  the  former  decreases  and  the  latter  increases. 

number  of  runs  were  made  in  real  argon,  using  the  thermodynamic  table 
for  triply-ionized  argon  described  in  Ref.  2.2  which  is  valid  for  densities 
above  10~7  g/cm3  and  internal  energies  up  to  2.5  *  10 12  erg/g.  Below  10"7  g/ca3 
the  program  uses  the  perfect  gas  law  for  argon. 

The  nozzle  geometry  was  the  conical  shape,  with  a  10°  half  angle  and  a 
throat  diameter  of  0.1016  cm. 

The  location  of  the  density  and  internal  energy  changes  was  taken  at  2  cm 
from  the  nozzle  throat.  For  a  stagnation  temperature  of  300  K,  which  yields  a 
limiting  velocity  of  5.59  x  104  cm/s,  2  cm  is  traversed  in  36  vs.  This  condition 
was  one  chosen  for  the  second  pulse  experiments,  and  was  used  in  the  calculations  to 
simulate  those  experiments. 

The  stagnation  conditions  used  were  10  atm  and  300  K.  The  runs  were  made 
in  the  blast  wave  mode  of  the  program.  A  given  amount  of  laser  energy  is  put  into  a 
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small  region  of  the  nozzle  near  the  throat  to  determine  the  Initial  thermodynamic 
conditions  there,  and  then  the  calculation  proceeds  to  resolve  the  resulting  discon¬ 
tinuity,  much  like  a  shock  tube  flow.  For  these  runs,  1  J  was  put  into  the  region 
extending  from  the  throat  to  0.1  cm  downstream.  This  is  a  volume  of  1.125  x  10”3 
cm3.  The  density  there  was  taken  to  be  uniform  at  the  value  in  the  perfect  gas  noz¬ 
zle  flow  at  0.1  cm,  which  is  3.66  x  10"3  g/cm3.  The  resulting  internal  energy  for 
1  J  is  2.429  x  lo12  erg/g,  leading  to  a  temperature  of  73,000  K  and  a  pressure  of 
2.0  x  108  dynes/cm2  (1950  atm). 

The  density  reduction  factor  was  taken  as  unity  (no  reduction),  100  and 
1000  in  three  runs,  to  compare  the  effects  of  having  the  tail  of  the  previous  pulse 
in  the  nozzle.  In  order  to  see  which  density  reduction  factor  is  most  realistic, 
the  density  at  2  an  from  the  throat  was  plotted  vs.  time  for  each  of  these  three 
runs,  as  shown  in  Fig.  2.1.  The  density  at  that  location  when  the  nozzle  is  filled 
with  gas  from  the  stagnation  chamber  is  8.48  x  1 0“5  g/cm3,  shown  by  the  line  at  the 
top.  The  open  squares  are  for  the  run  with  no  density  reduction,  the  open  circles 
for  a  factor  of  100  reduction  and  the  filled  circles  a  factor  of  1000  reduction. 
Figure  A  shows  that  there  is  not  much  difference  in  the  density  history  for  the  re¬ 
ductions  of  100  and  1000.  They  both  reach  8.5  x  10"7  g/cm3  at  about  21  jjs,  and  ex¬ 
trapolation  indicates  they  would  reach  8.5  x  10”8  g/cm3  at  about  40  us.  For  the  no¬ 
reduction  run,  the  densities  are  slightly  higher,  reaching  these  values  at  24  us  and 
55  us,  respectively.  At  36  us,  which  it  takes  to  traverse  2  cm  at  the  limiting 
velocity,  the  density  at  two  centimeters  is  nearer  to  the  factor  of  1000  reduction 
than  the  factor  of  100,  which  indicates  that  the  former  is  perhaps  a  more  realistic 
simulation . 

Of  course,  in  the  actual  case,  the  properties  of  the  gas  in  the  tail  have 
a  distribution  which  is  different  from  that  obtained  by  merely  changing  the  nozzle 
flow  state  by  some  factor.  The  present  calculations  should  be  looked  upon  as  a 
qualitative  indication  of  the  effect  of  having  the  previous  pulse  tail  in  the 
nozzle,  not  a  precise  numerical  simulation. 

In  order  to  evaluate  the  effect  of  the  density  reduction  we  will  compare 
the  results  of  the  computer  calculations  for  the  no-density-reduction  and  factor- 
of-1000  reduction  cases. 
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Figure  2.2  shows  uaVg  vs.  x  distributions  at  various  times,  where  x  is  cen¬ 
timeters  downstream  of  the  throat,  uavt,  is  measured  in  units  of  1.8625  *  10s  cm/s, 
and  the  times  of  the  various  distributions  are  in  us.  Figure  2.2a  is  for  no  density 
reduction,  while  Fig.  2.2b  is  for  the  factor  of  1000  reduction  at  2  cm.  The  reduced 
density  case  shows  considerably  higher  values.  To  convert  uaVg  to  specific  impulse 
in  seconds,  one  multiplies  uavg  by  190. 

Recall  from  Eq.  (2.1)  that  at  any  x  and  t,  the  value  of  uavg  is  the  ratio 
of  momentum  flux  passing  that  x,  integrated  up  to  t,  to  the  same  integral  for  the 
mass  flux.  Thus,  uaVg  depends  on  both  x  and  t.  Until  the  shock  arrives  at  any 
station,  the  density  is  fixed  and  the  velocity  is  approximately  the  limiting  veloc¬ 
ity,  which  is  5.5875  x  104  cm/s  for  300  K.  This  is,  in  fact,  the  value  that  uavg 
takes  on  in  front  of  the  shock,  as  Fig.  2.2a  shows.  When  the  shock  passes,  the  gas 
is  accelerated,  and  uaVg  increases. 

In  the  no-density-reduction  case,  the  shock  keeps  swallowing  significant 
amounts  of  mass  as  it  moves,  and  this  mass  must  be  accelerated.  In  fact,  the  mass 
per  unit  length  the  shock  swallows  is  nearly  constant,  since  it  is  the  mass  flow  in 
the  nozzle  divided  by  the  undisturbed  gas  velocity,  which  is  nearly  the  limiting 
speed  for  most  of  the  nozzle.  As  a  result,  the  average  velocity  decreases  with  x  at 
any  given  time,  and  also  decreases  with  time  at  a  fixed  x,  after  the  shock  passes. 
Figure  2.2a  shows  that  at  10  cm,  shortly  after  the  shock  passes  (32.6  us)  the  X8p  is 
190  s.  At  16  cm  and  61.6  us,  the  ISp  has  dropped  to  170  s. 

When  there  is  a  density  reduction,  the  shock  swallows  only  one-thousandth 
as  much  gas  per  unit  length,  after  the  reduction,  and  so  can  move  faster  and  accel¬ 
erate  the  gas  to  a  higher  speed.  Figure  2.2b  shows  that  at  10  cm  and  13.8  us,  the 
shock  has  passed  and  the  ISp  is  722  s.  At  17  cm  and  21.8  us,  the  ISp  is  down  to 
646  s. 

The  increase  in  ISp  with  density  reduction  is  brought  on  by  the  consid¬ 
erable  reduction  in  mass  fluence.  This  effect  can  be  seen  in  the  distributions 
of  mass  fluence  presented  in  Fig.  2.3.  The  curves  for  no  density  reduction  are 
shown  in  Fig.  2.3a,  where  the  top  of  the  ordinate  scale  is  1.5  x  10-1  in  units  of 
3.14  x  10“3  g.  At  10  cm  and  32.6  us,  1.82  x  10~4  g  has  passed,  while  at  16  cm 
and  61.6  us,  1.48  x  10~4  g  has  passed.  The  station  of  maximum  mass  fluence  at  a 
given  time  is  some  distance  behind  the  shock  wave.  For  a  factor  of  1000  density 
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Fig.  2.2a  Mass  averaged  velocity  distributions  at  various  times.  No  density 

reduction.  Real  argon  for  stagnation  conditions  of  300  K  and  10  atm 
with  1  J  deposited  between  the  throat  and  0.1  cm. 
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Fig.  2.2b  Mass  averaged  velocity  distributions  at  various  times.  Density  re¬ 
duced  by  a  factor  of  1000  downstream  of  2  cm.  Real  argon  for  stagna¬ 
tion  conditions  of  300  K  and  10  atm,  with  1  J  deposited  between  the 
throat  and  0.1  cm. 


14 


5527-4 


X(cm) 


Pig.  2.3a  Mass  fluence  distribution  at  various  times.  No  density  reduction. 

Real  argon  for  stagnation  conditions  of  300  K  and  10  atm,  with  1  J 
deposited  between  the  throat  and  0.1  cm. 
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Fig.  2.3b  Mass  fluence  distributions  at  various  times.  Density  reduced  by  1000 
at  downstream  of  2  cm.  Heal  argon  for  stagnation  conditions  of  300  K 
and  10  atm,  with  1  J  deposited  between  the  throat  and  0.1  cm. 
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reduction  the  corresponding  curves  are  shown  in  Pig.  2.3b  where  the  top  of  the  or¬ 
dinate  scale  is  3  x  10“2  in  units  of  3.14  x  10“3  g  (which  is  one-fifth  the  value 
in  Fig.  2.3a).  Here  the  shape  of  the  mass  fluence  curves  is  much  different  than 
in  Fig.  2.3a  reflecting  the  fact  that  nearly  all  the  mass  started  out  between  x-0 
and  x*2.  This  amount  of  mass,  8.98  x  10"5  g,  ha6  the  value  2.86  x  10“2  on  the 
ordinate  of  Fig.  2.3b,  and  the  figure  shows  that  at  23.8  us,  when  the  shock  is 
nearly  at  20  an,  only  90%  of  that  amount  of  mass,  at  most,  has  passed  any 
station.  Some  of  the  remaining  10%  is  still  in  the  initial  part  of  the  nozzle, 
while  some  may  represent  numerical  error,  since  the  time  integrations  at  each 
station  were  done  rather  crudely.  The  amounts  of  mass  which  have  passed  10  an 
at  13.8  us,  or  17  cm  at  21.8  us  are  lost  at  the  bottom  of  the  ordinate,  but  are 
less  than  3  x  10~®  g.  It  is  these  small  values  which  are  the  denominator  of 
uavg  at  the  corresponding  points  in  Fig.  2.2b. 

The  numerator  of  uaVg,  the  momentum  fluence,  is  presented  in  Fig.  2.4. 

For  no  density  reduction.  Fig.  2.4a,  the  top  of  the  ordinate  is  0.15  in  units  of 
5.85  x  102  dyne-s,  while  the  top  of  the  ordinate  in  Fig.  2.4b  is  0.1  in  the  same 
units.  Again,  the  values  for  the  no-density-reduction  case  are  higher  than 
those  for  the  factor-of-1000  reduction.  However,  this  reduction  is  not  as  large 
as  for  the  mass  fluence,  so  the  momentum  to  mass  fluence  ratio,  uavg,  is  higher 
in  the  density-reduced  case,  as  we  have  already  seen  in  Fig.  2.2. 

Figure  2.5  shows  that  the  density-reduced  case  accelerates  the  gas  to 
higher  speeds  than  the  no-reduction  case.  The  gas  speed  for  the  no-reduction  case 
is  presented  in  Fig.  2.5a  where  the  top  of  the  ordinate  scale  is  2  in  units  of 
1.8625  x  10^  cm/s.  The  rapid  oscillations  up  to  5  cm  go  with  the  61.6  and  62.9  us 
distributions.  They  are  caused  by  the  discontinuity  in  the  gas  law  at  a  density 
of  10~7  g/cm3  mentioned  above,  where  real  argon  is  used  above  10~7  and  perfect  ar¬ 
gon  below  10-7.  Unfortunately,  from  x*1  to  5  cm  the  density  is  very  close  to  10“7 
at  those  times,  which  produces  the  large  oscillations  in  the  velocity  as  well  as 
all  the  thermodynamic  variables.  However,  Fig.  2.5a  nevertheless  shows  gas  speeds 
after  the  shock  in  the  range  of  0.9  to  1.6.  The  corresponding  curves  for  the 
density-reduced  case  are  shown  in  Fig.  2.5b,  where  the  top  of  the  ordinate  scale 
is  5  in  units  of  1.8625  x  10s  cm/s  (which  is  2.5  times  the  value  in  Fig.  2.5a. 

The  gas  speeds  after  the  shock  are  in  the  range  of  3.3  to  4.5,  about  three  times 
as  high  as  for  the  no-reduction  case. 
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Fig.  2.4b  Momentum  fluence  distributions  at  various  times.  Density  reduced  by 
a  factor  of  1000  downstream  of  2  cm.  Real  argon  for  stagnation  con¬ 
ditions  of  300  K  and  10  atm.  with  1  J  deposited  between  the  throat 
and  0.1  cm. 
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Fig.  2.5a  Gas  velocity  distributions  at  various  times.  No  density  reduction. 

Real  argon  for  stagnation  conditions  of  300  K  and  10  atm,  with  1  J 
deposited  between  the  throat  and  0. 1  cm. 
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The  tick  marks  on  the  curves  in  Fig.  2.5b  give  the  location  at  which  the 
mass,  integrated  from  the  throat,  equals  the  mass  originally  located  from  x>0  to 
2  an,  namely  8.98  *  10-5  g.  This  is  the  gas  which  represents  the  new  gas  in¬ 
jected  from  the  stagnation  chamber  for  the  current  pulse.  The  gas  from  the  tick 
marks  to  the  corresponding  shock  fronts  was  the  reduced-mass  gas  originally 
downstream  of  2  cm,  meant  to  represent  the  gas  remaining  from  the  previous 
pulse.  The  figure  shows  that  both  the  residual  gas  and  the  newly  injected  gas 
have  been  accelerated  to  high  speeds. 

This  simulation  shows  that  if  we  fill  the  nozzle  only  partially  with 
fresh  gas,  and  fire  the  laser  while  the  tail  of  the  gas  from  the  previous  pulse 
is  still  in  the  nozzle,  we  can  impart  a  higher  momentum  per  unit  mass  to  the  gas, 
and  thus  a  higher  specific  impulse. 

As  an  aid  to  the  experimental  diagnostics,  plots  were  also  made  of  the 
pressure  vs.  time  at  two  stations  in  the  nozzle  where  pressure  transducers  were 
installed,  namely  4.4  cm  and  9.4  cm.  These  plots  are  shown  in  Fig.  2.6,  where 
Fig.  2.6a  is  for  the  no-reduction  case,  and  Fig.  2.6b  for  the  factor  of  1000 
density  reduction.  In  the  former  case,  we  see  that  the  pressure  jumps  to  its 
highest  value  when  the  shock  arrives,  and  then  decays,  for  both  stations.  For 
the  reduction  case.  Fig.  2.6b,  the  pressure  behavior  is  somewhat  different.  The 
pressure  at  4.4  cm  jumps  when  the  shock  arrives,  but  after  a  small  wiggle,  con¬ 
tinues  to  rise  for  some  time  before  decaying.  The  pressure  at  9.4  cm  also  jumps 
upon  shock  arrival,  then  decays,  but  reaches  a  minimum  and  rises  again.  (The 
dashed  part  of  the  curve  is  an  interpolation  necessitated  by  an  error  in  specify¬ 
ing  the  location  of  the  pressure  station  during  part  of  the  calculation,  so  the 
depth  of  the  reduction  is  not  known.  However,  the  right  end  is  defined  by  two 
correctly  calculated  points,  so  there  is  no  doubt  the  curve  rises.) 

The  behavior  of  the  pressure  shown  in  Fig.  2.6b  has  some  implications  for 
the  interpretation  of  the  pressure  transducer  data  obtained  in  the  double  pulse 
experiments  (see,  for  example.  Fig.  3.3).  In  those  experiments,  the  pressure 
peak  that  is  observed  at  each  transducer  following  the  second  laser  pulse  is 
taken  to  indicate  the  arrival  of  the  slug  of  propellant  gas  heated  by  that  pulse. 
The  propellant  exhaust  velocity  is  then  inferred  from  the  transit  time  of  this 
"wave"  between  the  two  pressure  stations.  It  appears  from  Fig.  2.6b  that  this 
interpretation  may  have  to  be  examined  more  closely. 
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Fig.  2.6b  Pressure  history  at  pressure  transducer  locations.  Density  reduced 
by  a  factor  of  1000  downstream  of  2  cm.  Real  argon  for  stagnation 
conditions  of  300  K  and  10  atm,  with  1  J  deposited  between  the  throat 
and  0.1  cm. 
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Calculations  for  Blast  Waves  in  Quiescent  Air 


As  an  aid  to  the  interpretation  of  some  experiments  in  quiescent  air, 
three  calculations  were  performed.  In  all  cases,  the  initial  state  of  the  gas 
in  the  nozzle  was  1  atm,  300  K,  with  zero  velocity.  Energy  was  deposited  into 
the  zone  near  the  conical  nozzle  throat.  The  gas  there  was  taken  to  be  at  the 
same  density  as  the  nozzle  gas,  1.1719  x  10”3  g/cm3,  but  at  4000  K.  The  size  of 
the  energy  deposition  zone  was  varied  to  allow  different  amounts  of  energy  to  be 
deposited.  In  one  run,  the  region  from  the  throat  to  3  cm  was  taken  as  hot, 
corresponding  to  a  deposition  of  4  J.  In  a  second  run,  the  region  to  1.8  cm  was 
hot,  corresponding  to  a  deposition  of  1  J.  In  both  cases,  perfect  air, 
was  used. 


The  information  desired  from  these  runs  was  the  trajectory  of  the  shock 
wave,  and  its  peak  pressure.  Correlation  formulas  for  these  quantities  were  ob¬ 
tained  in  a  manner  similar  to  that  described  in  Section  3.14  of  Ref.  2.2,  using 

the  appropriate  similarity  theory  of  blast  waves  into  quiescent  air.  The  length 
2  1/3 

scale  is  (E/p0a0),  where  E  is  the  energy  deposited,  and  pc  and  Bq  are  the  den¬ 
sity  and  sound  speed  in  the  quiescent  air.  (The  latter  two  for  the  present  case 
are  1.1719  x  10“3  g/cm3  and  3.479  x  lo4  cm/s.)  The  time  scale  is  this  length 
divided  by  Sq.  Therefore,  the  similarity  variables  are 

£  -  x/(E/p  a  V{3  T  -  t/{E/p  aV/3  (2.2) 

o  o  o  o 


for  length  and  time. 

The  shock  trajectory  in  these  variables  was  calculated  for  each  run. 

The  location  of  the  shock  at  each  time  was  defined  by  the  maximum  pressure  point, 
and  the  corresponding  value  of  5  was  denoted  as  £p.  A  plot  of  £p  vs.  t  for  the 

4  J  and  1  J  cases  are  shown  as  the  circles  and  triangles  in  Fig.  2.7.  It  can  be 

seen  that  the  scaling  represented  by  Eq.  (2.2)  reduces  the  two  trajectories  to 
the  same  curve,  in  spite  of  the  factor  of  4  difference  in  energy. 

In  the  experiments  to  be  described  later,  probes  were  at  x**4  and  9  cm, 

which  translate  to  a  pair  of  values  of  £  for  each  of  the  two  energies.  These 
values  are  denoted  by  horizontal  ticks  on  Fig.  2.7.  Only  the  values  of  £  between 
these  marks  are  interesting  in  using  the  data  from  the  two  probes. 
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To  correlate  these  points#  a  straight  line  in  the  log-log  plot  of 
Fig.  2.7  was  fitted  to  the  points  by  least  squares#  excluding  the  three  lowest 
circle  points  and  the  one  lowest  triangle.  This  portion  of  the  plot  is  in  the 
region  where  the  calculation  is  just  starting,  at  early  times#  and  shows  some 
deviation  from  a  straight  line.  The  equation  of  the  line  is 


£p  -  2.76  t°* 55926 


(2.3) 


When  we  revert  to  physical  variables, 
E  «  7.16  x  10-15  x 


using  Eq.  (2.2), 

6.807^3.807 

P 


and  solve  for  E,  we  find 
(2.4) 


in  cgs  units.  This  relation  expresses  the  energy  in  terms  of  the  time  of  shock 
arrival  at  the  station  Xp,  and  can  be  used  to  infer  the  energy  from  the  measure¬ 
ment  of  pressure  peak  arrival  time. 


A  similar  relation  between  the  peak  pressure  and  the  location  of  this 
peak  can  be  obtained.  For  the  pressure,  blast  wave  theory  suggests  that  the 
similarity  variable  is 


n  -  (p  -p  )/p  a*  (2.5) 

p  p  o  o  o 

where  pQ  is  the  pressure  in  the  quiescent  gas#  which  is  1  atm  in  the  present 
case.  A  plot  of  lip  vs.  £  for  the  4  and  1  J  runs  is  presented  as  the  circles  and 
triangles  in  Fig.  2.8.  Again  the  two  runs  fall  on  a  single  curve  except  at  low 
values  of  £p,  which  correspond  to  early  times.  The  vertical  ticks  again  are  the 
locations  of  the  two  probes  for  the  two  runs. 

A  straight  line  was  fitted  to  the  circles  and  triangles,  deleting  the 
three  circles  and  two  triangles  at  the  upper  end.  The  fit  is 


it  e  . a » — 1 .918 
II  «  6.6103  £ 

P  P 


which  translates  to 


E 


2.55  x  IQ4 


•(2.6) 


(2.7) 
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in  cgs  units.  This  relation  can  be  used  to  infer  the  energy  front  a  measurement 
of  peak  shock  pressure. 

A  third  run  into  quiescent  air  was  made  as  a  check  case,  to  compare  with 
a  run  made  by  Brode  in  the  early  days  of  the  use  of  computers  for  fluid  mechanics 
calculations,  Ref.  2.3.  He  calculated  the  blast  wave  resulting  from  the  explosion 
of  a  sphere  of  air  (Y“1*4)  at  ambient  density  but  a  pressure  of  121  atm,  into  an 
ambient  gas  at  1  atm.  He  used  Lagrangian  coordinates  and  a  calculation  algorithm 
introducing  artificial  viscosity,  in  contrast  to  the  Eulerian  coordinates  and  lax 
calculation  algorithm  used  here. 

This  case  has  an  energy  deposition  which  corresponds  to  approximately  9  J 
in  the  10°  nozzle  used  for  the  present  calculations.  The  results  of  this  calcula¬ 
tion,  in  similarity  variables,  are  shown  as  the  squares  on  Figs.  2.7  and  2.8.  On 
the  former,  the  squares  deviate  only  slightly  from  the  curve  defined  by  the  cir¬ 
cles  and  triangles.  In  Fig.  2.8,  the  deviation  is  larger.  However,  it  must  be 
remembered  that  the  initial  pressure  in  the  4  and  1  J  runs  was  only  13.3  atm,  com¬ 
pared  to  121  in  the  Brode  run.  It  also  appears  from  Fig.  2.8  that  the  squares  are 
merging  with  the  circles  and  triangles  as  the  pressure  drops. 

The  results  of  the  present  calculation  for  the  Brode  case  were  compared 
with  a  small-scale  plot  of  peak  pressure  vs.  distance  in  Ref.  2.3.  The  agreement 
seemed  to  be  satisfactory,  lending  additional  confidence  to  the  results  of  the 
present  computer  program. 
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3.  EXPERIMENTS 


This  section  describes  the  experimental  measurements  that  were  made  to 
assess  the  performance  of  a  pulsed  laser-heated  thruster.  For  the  most  part, 
the  experiments  that  were  carried  out  were  very  similar  to  those  performed 
previously3,1  -  both  in  measurement  technique  and  method  of  analysis.  The  ex¬ 
periments  were  performed  using  a  small  scale  conical  rocket  nozzle  and  either 
argon  or  hydrogen  propellant.  The  propellant  gases  were  heated  by  pulses  of  a 
TEA  CO2  laser  (X  *  10.6  ym,  Ep  »  10  J)  and/or  an  e-beam  pumped  XeF  laser  (X  * 
0.35  um,  Ep  *  4  J).  The  laser  radiation  was  focussed  into  the  nozzle  using 
external  focussing  mirrors  and  concentrated  to  sufficient  intensity  to  achieve 
propellant  breakdown  just  downstream  of  the  throat.  The  resulting  dynamics  of 
the  laser-generated  propellant  shock  waves  were  monitored  using  high  speed 
piezoelectric  pressure  transducers  mounted  in  the  nozzle  wall.  Using  the  nu¬ 
merical  code  simulations  described  in  the  previous  section,  the  pressure 
transducer  data  were  analyzed  to  infer  estimates  of  rocket  performance.  The 
performance  parameters  evaluated  include  conversion  efficiency  of  laser  energy 
to  blast  wave  energy,  laser  energy  to  exhaust  kinetic  energy,  and  specific  im¬ 
pulse,  i.e.,  mass  averaged  exhaust  velocity.  All  experiments  were  performed 
with  the  nozzle  in  a  vacuum  background  (P  <  10”4  atm) . 

The  following  subsections  describe  the  experimental  apparatus  and 
technique,  the  rocket  performance  data  and  its  analysis,  and,  finally,  the 
conclusions  that  can  be  drawn  from  the  experimental  results. 

3. 1  Experimental  Apparatus 

A  schematic  diagram  of  the  overall  experimental  apparatus  is  shown  in 
Fig.  3.1.  In  the  following  subsections,  each  major  component  of  the  apparatus 
will  be  described. 

3.1.1  Lasers 

The  Lumonics  K-103  pulsed  CO2  TEA  laser  was  Installed  in  a  carefully 
designed  and  electrically  shielded  metal  box  in  order  to  minimize  EMI  noise 
pickup  by  the  electronic  diagnostics.  Attention  was  paid  .<  shielding  against 
rf  radiation,  avoiding  ground  loop  problems,  and  filtering  the  power  line 
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inputs.  The  CO2  laser  was  placed  as  close  to  the  laser  propulsion  test  chamber 
as  space  allowed,  thus  minimizing  the  spreading  of  the  beam  as  it  travelled 
along  the  optical  train  and  obviating  the  need  for  any  beam  reducing  optics. 

The  details  of  the  pulsed  XeF  laser  (0.35  um)  have  been  described  pre¬ 
viously.^*^ 


3.1.2  Test  Chamber  and  Optical  Train 

The  vacuum  test  chamber  used  was  the  same  as  that  described  previous¬ 
ly. 3*  1  As  shown  in  Fig.  3.1,  a  KC1  window  (2"  diameter  and  1/8“  thickness) 
was  placed  on  the  plexi-glass  flange  to  allow  transmission  of  the  CO2  laser 
pulse  (10.6  ym)  into  the  vacuum  test  chamber.  The  CO2  laser  beam  was  then 
turned  90°  by  a  two-by-two  inch  gold  surface  mirror  and  directed  to  an  8" 
aluminum  coated  concave  surface  mirror  (f  *  160  cm).  The  two-by-two  inch 
folding  mirror  was  placed  in  such  a  way  that  it  fit  inside  the  hole  of  the 
square  annular  XeF  laser  beam. 

The  XeF  laser  was  directed,  as  shown  in  Fig.  3.1,  to  the  8"  mirror  de¬ 
scribed  in  the  above.  The  two  laser  beams  were  then  directed  onto  the  3" 
aluminum  coated  off-axis  parabola.  This  off-axis  parabola  mirror  focused  the 
two  laser  beams  approximately  21  cm  from  the  mirror.  The  focal  location  of 
the  CO 2  laser  pulse  was  carefully  identified  by  observing  breakdown  sparks  in 
air.  For  the  XeF  laser  a  He-Ne  laser  beam  precisely  adjusted  to  match  the 
profile  and  propagation  axis  of  the  XeF  beam,  was  used  to  determine  the  focal 
point  by  observing  the  He-Ne  image.  It  was  found  that  the  focal  plane  of  the 
CO2  laser  beam  was  approximately  2.5  cm  beyond  that  of  the  XeF  laser.  This 
was  the  result  of  the  finite  wavefront  curvature  of  the  CO2  beam.  Neverthe¬ 
less,  when  the  apex  of  the  rocket  cone  was  placed  0.6  cm  beyond  the  focal 
point  of  the  XeF  laser,  the  CO2  laser  was  powerful  enough  to  induce  breakdown. 
As  a  result,  this  configuration  was  used  throughout  the  experiments.  * 

To  assess  how  much  optical  energy  was  entering  in  the  rocket,  it  was 
necessary  to  evaluate  the  reflectance  and  transmittance  of  the  various  mirrors 
and  windows  in  the  optical  train.  For  the  optical  train  associated  with  the 
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XeF  laser  beam,  the  following  values  were  used:  quartz  window  -  0.92  trans¬ 
mittance  and  0.08  reflectance;  aluminum  coated  surface  mirrors  -  0.88  reflec¬ 
tance.  The  resulting  transmission  was  0.655.  For  the  COj  laser  optical 
train,  the  following  values  were  used:  gold  mirror,  R  *  0.97;  KC1  window, 

T  -  0.92;  aluminum  coated  surface  mirrors,  R  “  0.95.  Resulting  transmission 
was  0.78. 


3.1.3  Rocket  Nozzle  Design  and  Propellant  Feed  System 

A  schematic  diagram  of  the  rocket  exhaust  cone  and  the  nozzle  is  shown 
in  Fig.  3.2.  The  rocket  was  made  of  aluminum  and  the  inside  surfaces  were 
polished  using  conventional  polishing  and  buffing  techniques.  The  throat  di¬ 
ameter  was  1  mm.  Piezoelectric  pressure  transducers  mounted  flush  with  the 
inside  wall  surface  were  at  4.3  cm  and  at  9.3  cm  from  the  apex  of  the  cone, 
respectively. 

For  most  of  the  runs,  argon  gas  was  used  as  the  propellant  and  gas 
flow  was  controlled  by  a  fast  acting  electrically  activated  solenoid  valve. 

The  solenoid  valve  maintained  gas  flow  for  approximately  25  ms  after  it  was 
activated,  thus  minimizing  the  pressure  rise  in  the  time  vacuum  chamber.  A 
piezo  pressure  transducer  was  also  mounted  on  the  rocket  plenum  to  monitor  the 
pressure  history  as  well  as  the  establishment  of  steady  state  gas  flow  through 
the  rocket  nozzle.  This  piezo  pressure  transducer  was  calibrated  against  a 
conventional  gas  regulator  pressure  gauge.  The  calibration  was  accomplished 
by  installing  the  same  solenoid  valve  and  plenum  assembly  (without  nozzle) 
just  after  the  regulator,  and  then  pulsing  the  gas  into  the  plenum.  In  this 
way,  pressure  drop  between  the  regulator  and  plenum  was  minimal.  For  the  ac¬ 
tual  socket  experiments,  in  which  5  feet  of  6  mm  i.d.  tubing  was  required  be¬ 
tween  the  regulator  and  solenoid  valve,  the  maximum  plenum  pressure  was  found 
to  be  approximately  90%  of  the  regulator  delivery  pressure. 

3.2  Experimental  Diagnostics 

3.2.1  Piezo  Pressure  Transducer 

Two  Kistler  pressure  transducers  were  used  to  monitor  the  laser- 
induced  shock  wave  arrival  time  to  calculate  shock  wave  parameters  and  the 
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post  shock  pressure  history.  The  piezo  pressure  transducers  were  mounted  In 
Teflon®  plugs  to  minimize  mechanical  vibrations,  and  they  were  mounted  flush 
with  the  Inside  surface  of  the  conical  rocket  wall.  The  two  piezo  pressure 
transducers  were  Installed  at  4.3  cm  and  at  9.3  cm  from  the  apex  of  the  cone. 
The  response  time  was  nominally  1  Us  and  the  responslvity  was  25  mV/psl. 

3.2.2  Laser  Energy  Measurement 

The  XeF  laser  energy  output  was  monitored  by  a  Lumonlcs  500  calori¬ 
meter  measuring  the  reflected  beam  from  the  quartz  entrance  window  (total  re¬ 
flectivity  of  8%)  as  shown  in  Fig.  3.1. 

Energy  measurements  were  carried  out  always  simultaneously  with  the 
laser-induced  shock  transit  time  measurements,  since  the  energy  ouput  of  the 
XeF  laser  varied  from  run  to  ruin.  This  allowed  the  evaluation  of  blast  wave 
energy  or  kinetic  energy  conversion  efficiency  of  the  XeF  laser  energy. 

C02  TEA  laser  (10.6  um)  energy  output  was  monitored  using  the  same 
calorimeter,  in  a  separate  measurement  from  the  transit  time  measurement, 
since  the  C02  laser  energy  output  was  very  consistent  from  run  to  run. 

The  C02  laser  beam  was  directed  to  the  calorimeter  through  an  8.75% 
transmittance  beam  attenuator.  The  observed  energy  output  was  found  to  be 
consistently  11.9  J. 

3.3  Experimental  Technique  and  Procedure 

The  test  chamber  was  evacuated  to  approximately  2  x  10'"3  torr.  For 
most  of  the  experiments,  argon  gas  was  used  as  the  propellant  and  was  pulsed 
by  a  fast  solenoid  valve  to  establish  steady  state  flow  conditions  in  the 
rocket  nozzle.  17  ms  was  required  to  reach  the  steady  state  condition  over 
the  delivery  pressure  range  of  from  50  psig  to  150  psig.  The  C02  TEA  laser 
was  then  pulsed  to  create  a  breakdown  in  the  argon  gas  and  to  expand  the  gas 
outside  of  the  rocket  cone.  The  breakdown  near  the  apex  was  of  sufficient 
pressure  to  temporarily  stop  the  argon  gas  flow  through  the  nozzle.  After  a 
predetermined  time  delay,  the  XeF  laser  was  fired  into  the  rocket.  The  pre- 
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determined  time  between  laser  pulses  was  chosen  to  be  short  enough  so  that 
propellant  gas  refilled  only  a  small  portion  of  the  nozzle.  As  the  XeF  laser 
was  fired,  a  simultaneous  measurement  of  energy  was  carried  out  as  described 
in  the  previous  section.  Pressure  wave  transit  times,  and  pressure  history  in 
the  plenum  and  in  the  rocket  cone  were  recorded  by  photographing  oscilloscope 
traces  of  the  transducer  outputs. 

In  the  double  pulse  experiments,  the  purpose  of  the  first  laser  pulse 
is  simply  to  generate  a  nozzle  flow  field  for  the  second  pulse  that  simulates 
what  would  be  encountered  in  an  RP  thruster.  As  described  above,  the  break¬ 
down  due  to  the  first  pulse  creates  an  overpressure  downstream  of  the  throat 
which  temporarily  stops  the  flow  from  the  stagnation  chamber.  Then,  as  the 
blast  wave  expands  down  the  nozzle  and  the  pressure  decays  sufficiently,  the 
flow  restarts  and  the  nozzle  refills.  The  extent  to  which  the  nozzle  refills 
depends  on  the  delay  time  to  the  subsequent  laser  pulse.  Thus,  at  the  time 
the  second  (or  Nth)  laser  pulse  enters  the  rocket,  the  region  of  the  nozzle 
near  the  throat  contains  fresh,  cold  gas  which  has  entered  through  the  throat. 
However,  the  rest  of  the  nozzle  contains  heated,  low  density  gas  left  over 
from  the  previous  pulse.  It  is  this  flow  field  situation  which  is  modeled  in 
Sec.  2.  In  the  subsequent  data  analysis,  experimental  values  of  specific  im¬ 
pulse  and  laser  to  propellant  kinetic  energy  conversion  efficiencies  are  de¬ 
rived  from  second  pulse  pressure  transducer  data. 

Other  experiments  were  carried  out  to  determine  the  conversion  effici¬ 
ency  of  laser  energy  to  blast  wave  energy.  For  these,  only  one  laser  was 
fired,  either  the  XeF  or  the  COj  laser,  and  the  same  measurement  procedures 
were  performed. 

The  experimental  test  matrix  is  given  in  Table  3.1.  The  time  between 
laser  pulses  was  varied  from  30  to  100  us.  For  the  shortest  time  interval, 
the  shock  due  to  the  deposition  of  CO 2  laser  energy  was  still  in  the  nozzle. 
The  longest  time  delay  was  limited  by  signal  to  noise  and  the  lower  efficiency 
of  converting  deposited  energy  to  directed  kinetic  energy.  The  time  interval 
between  laser  pulses  was  controlled  by  a  California  Avionics  digital  time  de¬ 
lay  generator. 
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TABLE  3.1 


Test  Matrix  of  Laser  Propulsion  Experiments 


Footnote:  Numbers  in  the  test  matrix  indicate  mass  of  propellant  (in  grams) 

heated  by  the  XeF  laser  pulse.  The  propellant  used  for  the  above  ex¬ 
periment  was  argon. 
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4  Experimental  Results  and  Data  Reduction 
3.4.1  Data  Reduction  and  Error  Analysis 

Representative  oscilloscope  records  for  the  measurements  of  the  XeF 
laser  energy,  pressure  in  the  plenum,  and  shocks  in  the  nozzle  are  shown  in 
Figs.  3.3a,  b,  and  c. 

Peak-voltage  shown  in  Fig.  3.3a  is  proportional  to  the  XeF  energy  out¬ 
put  reflected  from  the  quartz  window.  The  sensitivity  of  Lumonics  50D  calo¬ 
rimeter  is  0.39  volt/J  and  20%  of  the  XeF  laser  beam  fell  outside  of  the  calo¬ 
rimeter,  thus  the  XeF  energy  output  is  given  by 

E(J)  -  12.5  x  x  1.25  ,  (3.4.1) 

0  •  39 

where  V  is  the  observed  peak-voltage  shown  in  Fig.  3.3. 

Pressure  history  in  the  plenum  shown  in  Fig.  3.3b  indicates  that  after 
the  solenoid  valve  opened,  the  pressure  reaches  steady  state;  a  large  spike 
which  appeared  at  17  ms  shows  the  time  that  the  CO2  and  the  XeF  lasers  were 
fired.  The  steady  state  pressure  in  the  plenum  determines  the  argon  mass  flow 
rate  through  the  1  mm  diameter  throat,  which  is  given  by 


*  *  * 

p  u  A 


(3.4.2) 


-  P(atm) (0.2387) 

where  p*  is  the  propellant  density  in  the  throat,  u*  is  the  sonic  velocity  of 
argon,  and  A*  is  the  nozzle  throat  area.  Mass  introduced  into  the  nozzle  be¬ 
tween  the  laser  pulses,  dm  is  given  by 


Am  ■  m  *  At  ,  (3.4.3) 

where  At  is  the  time  between  the  two  laser  pulses.  Equation  (3.4.3)  'is  an 
upper  limit  to  the  expelled  mass  since  it  assumes  that  the  propellant  flow 
restarts  to  the  steady  state  value  immediately  after  the  C02  laser  was  fired. 
Calculations  carried  out  using  the  flow  simulation  code  described  in  Sec.  2.2 
suggest  that,  in  fact,  for  the  conditions  of  the  experiment,  flow  from  the 
plenum  may  not  restart  for  approximately  3  microseconds.  Thus,  for  laser  in¬ 
terpulse  times  ■  30  microseconds,  the  error  inherent  in  using  Eq.  (3.4.3)  is 
less  than  10%. 
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Fig.  3.3  Represer  *- utive  oscilloscope 
records  of  the  XeF  laser  energy 
measurement  (-a) ,  the  pressure 
history  in  the  plenum  (-b) ,  and 
the  laser-induced  shock  wave  meas¬ 
urement  (-c) .  (a)  50  ms/div,  50 

mV/div:  (b)  10  ms/div  100  mV/div; 

discontinuity  on  the  plateau  re¬ 
gion  indicates  the  time  when  the 
CC>2  and  the  XeF  lasers  were  fired: 
(c)  10  ps/div,  upper  500  mV,  lower 
100  mV/div,  respectively;  +  indi¬ 
cates  the  CO2  laser  firing;  *  and 
indicate  the  arrival  of  the  CO2 
laser-induced  blast  wave  at  the 
first  and  second  pressure  trans¬ 
ducer  station,  respectively;  'ff' 
indicates  the  XeF  laser  firing; 

A  and  5  indicate  the  arrival  of 
the  XeF  laser- induced  shock  wave 
at  the  first  and  second  pressure 
transducer  station,  respectively. 
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In  Fig.  3.3c,  the  following  events  were  identified  as  shown  with  signs: 
(1)  the  time  the  CO2  laser  was  fired  {  +  );  (2)  the  arrival  of  the  CO2  laser- 
induced  blast  wave  at  the  first  station  ( * )  and  at  the  second  station  ( * ) » 

(3)  the  time  the  XeF  laser  was  fired  (30  us  after  (DM'ft);  (4)  the  arrival  of 
the  XeF  laser-induced  shock  wave  at  the  first  (A)  and  second  stations  (&). 

The  CO2  (or  XeF)  laser  energy  deposited  in  the  steady  state  flow  of 
argon  is  given  by 


_  0.4142  *  P(atm)  *  {X(cm)}4'38  ..  .  .. 

E(J)  =  - — — -  (3.4.4) 

(t(us)}3,38 

where  P(atm)  is  the  pressure  in  the  plenum,  x(cm)  is  the  distance  between  the 
breakdown  point  and  the  piezo  pressure  transducer  stations,  and  t(us)  is  the 
transit  time  from  the  end  of  the  breakdown  to  the  arrival  of  the  blast  wave  at 
each  pressure  transducer  station.  The  derivation  of  formula  3.4.4  is  in  Ref. 
3.1.  Also,  for  experiments  carried  out  in  stationary  air,  the  laser  energy  de¬ 
posited  in  the  rocket  cone  is  given  in  terms  of  pressure  by 

E(erg)  -  2.41  x  10"3  »  pi. 5144  at  tjje  first  piezo  pressure  (3.4.5) 

transducer  station 


or 

E(erg)  «  2.83  x  10-2  x  pl.5144  at  the  second  piezo  pressure  (3.4.6) 

transducer  station 

where  P  is  the  observed  pressure  jump  from  the  pre-shock  condition  (in  dyne 
cm-2).  The  derivation  of  these  equations  is  given  in  Sec.  2.2. 

The  laser-induced  propellant  exhaust  velocity,  averaged  over  the  two 
piezo  pressure  transducer  stations,  w$s  obtained  by  measuring  a  transit  time 
between  the  two  pressure  spikes  as  shown  with  A  and  &  in  Fig.  3.3c.  This  anal¬ 
ysis  was  only  applied  to  the  second  laser  heating  pulse. 

Similarly  the  exhaust  kinetic  energy  of  the  propellant  was  taken  to  be 
K.E.  -  j  Am  •  v2  (3.4.7) 
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where  v  is  the  average  propellant  exhaust  velocity  (as  described  above)  and 
Am  is  the  mass  expelled  between  the  first  and  second  laser  pulses.  Specific 
impulse  is  taken  to  be  V/g  where  g  is  the  gravitational  constant. 

Since  the  events  of  interest  in  the  rocket  cone  take  place  on  a  time 
scale  comparable  to  the  response  time  of  the  piezo  pressure  transducers,  and 
over  distances  comparable  to  the  transducer  diameters  and  uncertainties  in  the 
breakdown  location,  it  is  important  to  estimate  the  maximum  errors  incurred  in 
the  experiments.  A  sample  error  analysis  of  the  data  in  Fig.  3.3c  is  given  be¬ 
low. 

Using  Eq.  (3.4.4)  the  maximum  error  in  determining  the  CO2  laser- 
induced  blast  wave  energy  is  given  by 


0.4142  *  {■ 


4.38  Px' 

3.38 


3.38  Px 


*  0.4142  x  {2.12  +  3.17  +  0.74}  (3.4.8) 

where  P  =  10.1  atm,  t  *  16  us,  and  x  *  9.3  cm  and  Ax,  At,  and  Ap  are  taken  to 
be  0.3  cm,  1.0  us,  and  0.5  atm.  Then,  the  absolute  value  of  AE  is  given  by 
2.5  J,  which  corresponds  to  1  i  40%  error.  As  seen  in  Eq.  (3.4.8),  the  ac¬ 
curacy  of  measuring  the  laser  energy  deposited  in  the  propellant  was  limited 
primarily  by  the  response  time  of  the  piezo  pressure  transducer  and  the  uncer¬ 
tainty  of  the  precise  breakdown  location  in  the  rocket  cone.  Similarly  the 
maximum  error  in  the  determination  of  the  propellant  kinetic  energy  resulting 
from  the  second  laser  pulse  is  given  by 
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where  x  »  5.08  cm,  P  *  10.1  atm,  t  “  10.5  us,  and  *  30  us  and  Atde^*y 

«  1  us.  Ax  ■  0.3  cm,  AP  «  0.5  atm,  and  At  -  1.5  us.  Then,  the  absolute  value 
of  AE  is  given  by  0.43  <7,  which  corresponds  to  -  50%.  As  seen  in  Eq.  (3.4.9), 
the  biggest  errors  arise  from  the  finite  dimension  of  the  pieso  pressure  trans¬ 
ducer  head  and  uncertainties  in  the  transit  time  measurements.  As  is  illus¬ 
trated  above,  the  systematic  errors  associated  with  the  present  measurements 
are  significant,  i.e.,  -  40%  for  the  blast  wave  energy,  ±  50%  for  the  exhaust 
kinetic  energy,  and  ±  20%  for  the  specific  impulse  (average  propellant  exhaust 
velocity).  To  reduce  these  errors,  improved  diagnostic  techniques  are  needed 
that  have  improved  spatial  and  temporal  resolution. 

3.4.2  Data  for  Energy  Conversion  of  Laser  Energy  to 
Blast  Wave  Energy 

Typical  oscilloscope  records  of  the  blast  wave  measurement  in  station¬ 
ary  1  atm  air  using  the  CO 2  and  the  XeF  lasers  are  shown  in  Figs.  3.4a  and  b, 
respectively.  For  these  experiments,  the  conversion  of  laser  energy  to  blast 
wave  energy  was  obtained  by  measuring  the  elapsed  time  from  breakdown  to  the 
arrival  of  the  blast  wave  at  each  station  and  measuring  the  pressure  jump  at 
each  station.  Equations  (3.4.5),  and  (3.4.6)  were  used  for  the  derivation  of 
the  blast  wave  energy.  The  experimental  results  are  shown  in  Table  3.2. 

Figures  3.5a  and  b  show  typical  oscilloscope  traces  used  to  determine 
the  conversion  efficiency  of  XeF  laser  energy  into  blast  wave  energy  in  an  ar¬ 
gon  flow.  In  these  experiments,  only  the  transit  time  from  breakdown  to  the 
arrival  of  the  blast  wave  at  the  second  station  was  employed  for  the  calcula¬ 
tion  of  the  conversion  efficiency.  The  conversion  efficiencies  achieved  using 
the  C02  and  the  XeF  lasers  are  plotted  as  a  function  of  mass  flow  rate  in 
Figs.  3.6  and  3.7.  It  is  seen  that  for  the  range  of  mass  flow  rates  investi¬ 
gated  the  coupling  of  the  C02  laser  (10.6  vm)  beam  with  argon  gas  is. relative¬ 
ly  Independent  of  the  flow  rate  while  the  coupling  of  the  XeF  laser  (0.35  inn) 
beam  appears  to  Increase  with  flow  rate.  Nevertheless,  for  mass  flow  rates  “ 
1.0  g/s,  both  lasers  appear  to  be  converted  to  propellant  blast  wave  energy 
with  an  efficiency  "  50  ±  20%. 
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TABLE  3.2 


Comparison  of  Two  Different  Types  of  Blast  Wave  Energy  Conversion 
Measurements  Using  the  CO 2  and  the  XeF  Lasers  (in  Stationary  Air) 


x  =  4.3  cm 
x  *  9.3  cm 


Observed  Energy  (J)  and  Conversion 
From  Transit  Time 

4.41  J  48% 

3.42  J  37% 


Efficiency  (%)  of  the  CO?  Laser 
From  Pressure  Jump 
6.49  J  70% 

5.19  J  56% 


x  *  4. 3  cm 
x  =  9.3  cm 


Observed  Energy  (J)  and  Conversion  Efficiency  (%)  of  the  XeF  Laser 
From  Transit  Time  From  Pressure  Jump 

0.677  J  38%  0.974  J  55% 

0.832  J  47%  1.27  J  72% 
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3.4.3  Data  for  Conversion  Efficiency  of  Laser  Energy  to  Gas 
Kinetic  Energy  Using  Two  Laser  Pulses  -  Determination 
of  Specific  Impulse 

Typical  oscilloscope  records  of  the  measurements  used  to  determine 
propellant  kinetic  energy  were  shown  In  Pigs.  3.3a,  b  and  c.  The  propellant 
kinetic  energy  appearing  in  the  exhaust  is  calculated  using  Eq.  3.4.7  with  the 
propellant  exhaust  velocity  determined  from  transit  time  measurements  of  the 
propellant  pressure  wave  induced  by  the  second  laser  pulse  (XeF  laser).  The 
reduced  data  are  plotted  in  Figs.  3.8-3.12  as  the  conversion  efficiency  of 

laser  energy  to  propellant  kinetic  energy  vs.  propellant  mass.  The  propellant 

mass  plotted  is  the  mass  expelled  between  the  first  and  second  laser  pulses, 

i.e.,  m  A  t,  where  At  is  the  laser  interpulse  time.  For  these  experiments, 

the  XeF  laser  energy  entering  the  rocket  was  typically  2.6  J. 

In  Fig.  3.13,  we  plot  the  average  conversion  efficiency  achieved  vs. 
argon  mass  flow  rate.  For  each  mass  flow  rate,  the  conversion  efficiency 
plotted  represents  an  average  for  the  various  laser  interpulse  times  studied. 
As  can  be  seen,  for  mass  flow  rates  £  1.6  g/s,  the  conversion  efficiency  is 
relatively  independent  of  mass  flow  rate  and  appears  to  take  on  an  average 
value  of  25  to  30%. 

From  Figs.  3.8  through  3.12  it  can  be  seen  that  specific  impulse  of  up 
to  500  s  was  achieved  in  argon  propellants  heated  by  the  pulsed  XeF  (0.35  um) 
laser.  Furthermore,  this  specific  impulse  was  obtained  with  conversion 
efficiencies  of  laser  energy  to  propellant  kinetic  energy  of  up  to  30  to  50%. 
These  results  compare  favorably  with  similar  results  that  were  presented  in 
Ref.  3.1  for  argon  propellants  processed  by  two  successive  10.6  pm  laser 
pulses.  In  the  latter  studies,  maximum  specific  impulses  of  500  s  in  argon 
and  1000  s  in  hydrogen  were  achieved  with  conversion  efficiencies  of  laser 
energy  to  exhaust  kinetic  energy  of  up  to  40  to  50%. 

Since  in  both  the  previous  studies3*1  and  the  present  work,  the  maxi¬ 
mum  specific  impulse  achieved  in  argon  was  500  seconds,  the  question  arose  as 
to  whether  this  represents  a  fundamental  physical  limit  in  argon.  However, 


49 


5514-8 


t 


in 

£  ro 


C  0 

O  O  Ml 

Cn  O  Vj 

H  ^  0) 

(0 


o>  —  e 

0) 

<0  O  0) 

e  o  j: 

m  -P 

0) 

D  n  o> 
w  Q) 
u  -P 
Q)  (/)  <0 
>  U 
CVH 
>.  01*0 
U  H  C 
C  *H 
<D  0 
•H  4J  M 
U  « 
•h  *C 
U-I  C 


U)  M 
COM 
C  M  W 
•H  Ml 
W  0 
U  U  • 


C  0)  vO 

o  c  • 

U  *H  fN 


CT'O  o 
V4  0) 

JC  >i 

C  to  O' 
Q)  m  U 
Q  0) 
u  c 

■H  01 
■P  • 

0)  01  Em 
C  \ 

•h  trx 

m  43 
01  CM  4J 


o> 

•H 

&4 


50 


ment,  which  was  evaluated  using  Eq.  (3.4.9) 
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of  standard  deviation. 
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because  the  earlier  experiments  were  performed  with  a  nozzle  geometry  signifi¬ 
cantly  different  from  the  present  one,  i.e.,  a  rapidly  expanding  paraboloid 
vs.  the  present  10*  cone,  it  was  difficult  to  draw  any  strong  conclusions. 
Thus,  it  was  decided  to  perform  two  pulse  experiments  in  the  conical  nozzle  in 
which  both  the  first  and  second  pulses  were  CO2  laser  pulses.  In  principle 
this  would  allow  direct  comparison  of  0.35  pm  and  10.6  pm  performance  under 
the  same  nozzle  conditions,  as  well  as  comparisons  of  the  performance  of  the 
parabolic  and  conical  nozzles  under  the  same  conditions. 

In  Fig.  3.14  we  plot  the  specific  impulse  data  obtained  for  two  10.6 
inn  laser  pulses  heating  argon  propellant  in  the  conical  nozzle.  As  was  de¬ 
scribed  above,  the  specific  impulse  is  determined  from  transit  time  measure¬ 
ments  of  the  nozzle  pressure  waves  induced  by  the  second  laser  pulse.  The 
data  of  Fig«  3.14  are  plotted  as  the  specific  impulse  observed  vs.  Am,  the 
propellant  mass  heated  by  the  second  CO2  laser  pulse.  As  before,  we  have 
taken  Am  to  be  equal  to  mAt,  where  At  is  the  time  delay  used  between  laser 
pulses.  The  diagonal  lines  plotted  with  the  data  represent  loci  for  different 
values  of  exhaust  kinetic  energy,  i.e.,  “  1/2  AmV2  *  1/2  Am(gISp)2.  As  can 

be  seen,  the  maximum  specific  impulse  achieved  in  these  experiments  was  ap¬ 
proximately  700  s.  For  a  propellant  mass  of  *  5  x  10”5  grams,  this  value  of 

specific  impulse  indicates  a  maximum  exhaust  kinetic  energy  of  ■  1.2  J.  With 

the  total  CO2  laser  energy  entering  the  rocket  being  approximately  6.5  J,  this 
corresponds  to  a  kinetic  energy  conversion  efficiency  of  18.5%.  It  is  inter¬ 
esting  to  note  that  a  specific  impulse  of  700  s  in  argon  is  40%  higher  than 
any  of  the  previous  results.  Comparing  conical  nozzle  results  only,  this  im¬ 
proved  performance  is  probably  due  to  the  higher  laser  energy  available  at 
10.6  pm  vs.  0.35  pm,  i.e.,  6.5  J  vs.  2.7  J.  Presumably  similar  results  could 

be  obtained  with  an  XeF  laser  if  the  laser  energy  were  increased. 

Based  upon  the  most  simple  arguments,  the  specific  impulse  is  expected 
to  scale  as  ~  ^Eiaser*  We  see  that  this  is  not  quite  the  case  when  we  compare 
the  specific  impulse  obtained  using  the  CO2  laser  and  that  achieved  with  the 
XeF  laser.  From  the  ratio  of  the  available  laser  energies,  ~  6. 5/2. 7,  the 
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PROPELLANT  MASS  (Am)  -  grams 

Fig.  3.14  Specific  impulse  data  obtained  for  argon  propellant  heated  by  two  C02 
TEA  laser  pulses.  Specific  impulse  is  determined  from  transit  time 
measurements  of  pressure  wave  induced  by  2nd  laser  pulse.  Dashed 
lines  plotted  with  data  indicate  loci  of  constant  propellant  kinetic 
energy  (>j  Am(g  Isp)2).  is  the  mass  expelled  through  throat 

between  first  and  second  laser  pulses. 


above  argument  would  lead  us  to  predict  a  specific  impulse  improvement  of  55%. 
The  observed  improvement  was  only  40%.  It  is  clear  that  in  future  studies  it 
would  be  worthwhile  to  investigate  how  the  achievable  specific  impulse  varies 
with  input  laser  energy. 

It  is  also  interesting  to  note  the  higher  specific  impulse  achieved 
with  the  present  conical  nozzle  compared  to  the  previous  parabolic  nozzle.  He 
believe  the  improved  performance  results  from  conditions  in  the  conical  nozzle 
being  more  conducive  for  efficient  laser  energy  deposition.  Apparently,  at 
the  lower  mass  flow  rates,  the  very  rapid  initial  divergence  of  the  parabolic 
nozzle  leads  to  downstream  gas  densities  that  are  too  low  to  support  strong 
LSD  wave  absorption . 

Finally,  measurements  similar  to  those  presented  for  argon  in  Fig. 

3.14  were  also  obtained  for  hydrogen  propellants.  The  results  are  given  in 
Fig.  3.15.  Note  the  significant  increase  in  specific  impulse  that  occurs  when 
hydrogen  is  used,  i.e.,  up  to  3000  seconds  in  hydrogen  vs.  a  maximum  of  700 
seconds  in  argon.  Qualitatively,  this  result  is  not  surprising  in  light  of 
the  much  lower  molecular  weight  of  hydrogen.  In  fact,  for  specific  impulse 
scaling  as  the  inverse  square  root  of  molecular  weight  (assuming  constant 
stagnation  temperature),  an  increase  of  /40/2  •  4.5  is  expected.  The  observed 
improvement  is  a  factor  of  ~  4.3. 

Comparing  the  results  in  Figs.  3.14  and  3.15  we  also  note  that  the 
conversion  efficiency  of  laser  energy  to  propellant  kinetic  energy  is  consid¬ 
erably  higher  in  hydrogen  than  in  argon.  In  fact,  at  maximum  specific  im¬ 
pulse,  the  conversion  efficiency  achieved  in  hydrogen  appears  to  be  approxi¬ 
mately  four  times  greater  than  in  argon,  i.e.,  ~  74%  vs.  ~  19%  in  argon. 
Because  the  data  in  Fig.  3.6  indicate  CO2  laser  energy  deposition  efficiencies 
in  argon  of  40  to  50%,  it  would  appear  that  the  low  overall  conversion  to  ex¬ 
haust  kinetic  energy  in  this  propellant  is  limited  primarily  by  loss  mech¬ 
anisms  other  than  imperfect  absorption.  Energy  loss  channels  that  are  likely 
to  be  important  are  plasma  radiation  losses  and  losses  to  residual  enthalpy  of 
the  exhaust  gas. 
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Fig.  3.15  Specific  impulse  data  obtained  for  hydrogen  propellant  heated  by  two  CO2 
TEA  laser  pulses.  Specific  impulse  is  determined  from  transit  time  meas 
sureroents  of  pressure  wave  induced  by  second  laser  pulse.  Dashed  lines 
are  loci  of  constant  propellant  kinetic  energy  ih  Am(g  ISp)^)*  Am  is 
the  mass  expelled  through  throat  between  first  and  second  la^ser  pifl?es. 


As  was  pointed  out  with  the  argon  results,  the  present  results  for  hy¬ 
drogen  using  the  conical  thruster  exhibit  a  significant  improvement  in  perfor¬ 
mance  over  the  previous  results  obtained  with  the  parabolic  nozzle . 3 » 1  In  the 
latter  experiments,  the  best  performance  obtained  with  hydrogen  propellant  was 
a  specific  impulse  of  1000  seconds  at  a  conversion  efficiency  of  laser  to  ex¬ 
haust  kinetic  energy  of  "  40%.  We  believe  the  better  performance  achieved  in 
the  present  investigations,  i.e.,  3000  seconds  with  conversion  efficiencies  up 
to  "  70%,  is  primarily  the  result  of  more  efficient  laser  energy  absorption  in 
the  conical  thurster.  In  the  previous  parabolic  thruster,  the  initial  nozzle 
divergence  was  so  large  that  attempts  to  increase  specific  impulse  by  reducing 
the  amount  of  propellant  mass  heated  led  to  nozzle  gas  densities  that  were  too 
low  to  support  efficient  laser  absorption. 


3.5  Summary  of  Results  and  Conclusions 

The  results  of  the  present  thruster  performance  experiments  can  be 
summarized  as  follows: 

(1)  Using  a  conical  nozzle  of  10*  half  angle  and  external  focusing 
optics  the  following  rocket  performance  parameters  have  been  achieved: 

•  Laser  to  blast  wave  energy  conversion  efficiencies  in  argon  propellant 
of  ~  40  to  60%.  Such  energy  deposition  efficiencies  were  achieved 
using  either  a  CO2  TEA  laser  (1  ■  10.6  urn,  Tp  ■  3  us,  Ep  ■  6.5  J)  or  a 
xenon  fluoride  excimer  laser  (X  *  0.35  um,  Tp  *  0.6  us,  Ep  *  2.7  J). 

•  Specific  impulse  in  argon  of  500  s  using  the  XeF  laser  (Ep  *  2.7  J) 
and  700  s  using  the  CO2  laser  (Ep  "  6.5  J).  These  specific  impulses 
were  achieved  with  corresponding  conversion  efficiences  of  laser  ener¬ 
gy  to  exhaust  kinetic  energy  of  up  to  40%  i  10%  and  19%  ±  5%,  respec¬ 
tively. 

•  A  maximum  specific  impulse  in  hydrogen  propellant  (using  the  CO2  laser 
as  the  heating  source)  of  ~  3000  seconds.  This  specific  impulse  was 
achieved  with  a  conversion  efficiency  of  laser  to  exhaust  kinetic 
energy  of  up  to  ~  70%. 

For  all  the  above  experiments,  specific  impulse  was  determined  from  transit 
time  measurements  of  the  pressure  wave  induced  by  the  second  laser  pulse  in  a 
two  pulse  heating  sequence.  Typical  interpulse  times  were  20  to  100  us. 
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(2)  Performance  obtained  with  the  conical  nozzle  and  CO2  TEA  laser 
pulses  represents  a  significant  improvement  over  similar  experiments  performed 
previously^* 1  using  a  parabolic  nozzle.  Previously  the  maximum  specific  im¬ 
pulse  that  could  be  obtained  was  500  s  in  argon  and  1000  s  in  hydrogen.  This 
is  to  be  compared  to  the  present  results  of  700  s  and  3000  s,  respectively. 

The  improvement  is  believed  to  result  from  better  laser  energy  absorption  that 
occurs  in  the  conical  nozzle  when  operating  at  low  mass  flow  rates. 

Based  upon  the  experimental  results  described  above,  the  following 
conclusions  can  be  made: 

( 1 )  There  appears  to  be  no  fundamental  reason  why  pulsed  laser-heated 
thrusters  can  not  be  operated  efficiently  at  visible/UV  laser  wavelengths. 

The  present  results  indicate  conversion  efficiencies  of  laser  energy  to  pro¬ 
pellant  energy  at  A  *  0.35  um  that  are  comparable  to  those  obtained  for  X  m 
10.6  pm.  The  only  significant  question  that  remains  with  respect  to  wave¬ 
length  scaling  is  whether,  if  a  self-focusing  nozzle  is  to  be  used,  sufficient 
focusing  can  be  achieved  at  the  short  wavelengths  to  achieve  rapid  laser- 
induced  breakdown.  Optical  design  considerations  for  a  self-focusing  nozzle 
are  discussed  in  Sec.  4. 

(2)  It  appears  that  the  values  of  specific  impulse  that  were  obtained 
previously  using  a  parabolic  nozzle, 3* 1  i.e.,  500  seconds  in  argon  and  1000 
seconds  in  hydrogen,  were  limited  by  poor  laser  energy  deposition  at  low  pro¬ 
pellant  masses.  Apparently,  the  very  rapid  initial  divergence  of  that  nozzle 
created  downstream  gas  densities  that  were  too  low  to  efficiently  absorb  the 
incident  laser  radiation.  In  contrast,  the  present  conical  nozzle,  with  its 
10*  half  angle  divergence,  eliminates  this  problem  by  maintaining  higher  gas 
density  downstream  of  the  throat.  The  most  dramatic  improvement  is  seen  in 
the  case  of  hydrogen  for  which  the  achievable  specific  impulse  increased  from 
1000  seconds  to  some  3000  seconds.  To  achieve  similarly  high  specific  inqmlee 
in  the  parabolic  nozzle  will  probably  require  operating  at  higher  mass  flow 
rates  and  with  significantly  larger  laser  pulse  energies.  Thus,  for  full 
scale  operation,  parabolic  nozzles  will  probably  not  present  a  limitation  to 
rocket  performance. 


(3)  Finally,  it  is  worth  noting  the  similarity  between  the  rocket  per¬ 
formance  parameters  we  have  obtained  for  a  pulsed  laser-heated  thruster  and 
those  typical  of  an  arc-heated  rocket  engine.  According  to  Refs.  3.2  and  3.3, 
arc-heated  rocket  engines  typically  operate  with  a  specific  impulse  in  the 
range  of  400  to  2000  seconds  and  a  maximum  propellant  temperature  of  ~  5500*C. 
The  typical  working  fluid  is  hydrogen.  The  internal  efficiency  of  such  en¬ 
gines,  i.e.,  fraction  of  electrical  input  energy  converted  to  kinetic  energy 
of  ejected  mass,  is  typically  15  to  20  percent,  (but  can  probably  be  increased 
to  some  50  to  60  percent  if  the  propellant  is  used  as  a  regenerative  coolant). 
The  similar  performance  observed  for  a  pulsed  laser-heated  thruster  is  not 
surprising  in  light  of  the  similarity  of  the  propellant  heating  process,  i.e., 
gas  breakdown  fund  plasma  formation  followed  by  a  thermodynamic  nozzle  ex¬ 
pansion. 


3.6  Recommendations  for  Future  Studies 

Perhaps  the  most  significant  result  of  the  present  studies  is  the  ex¬ 
perimental  evidence  for  achieving  a  specific  impulse  of  some  2500  to  3000  sec¬ 
onds  with  hydrogen  propellant.  Unfortunately,  because  this  result  depends  on 
a  theoretical  interpretation  of  the  laser-induced  pressure  waves  measured  in 
the  nozzle,  its  value  is  in  question.  For  verification,  additional  experi¬ 
ments  are  required  in  which  a  more  direct,  reliable  measurement  of  specific 
impulse  is  carried  out.  Ideally,  a  diagnostic  technique  is  needed  that  can 
provide  a  measurement  of  exhaust  gas  velocity. 

One  approach  which  appears  feasible  (particularly  for  exhaust  gas  vel¬ 
ocities  in  excess  of  106  cm/s,  i.e.,  ISp  >  1000  seconds)  is  to  measure  the 
Doppler  shift  of  spectral  line  emissions  emanating  from  the  exhaust.  In  hy¬ 
drogen,  a  good  spectral  line  to  monitor  might  be  the  well-known  red  line,  Ha, 
of  the  atomic  hydrogen  Balmer  spectrum.  For  particle  velocities  of  3  *  1 0 '• 
cm/s  d8p  "  3000  seconds),  the  Doppler  shift  would  be  vp/c  or  1  part  in  104. 
High  resolution  spectrometers  are  readily  available  which  can  resolve  1  part 
in  105. 


More  reliable  measurements  of  rocket  performance  would  also  be  made 
possible  by  performing  longer  duration  laser  tests.  Sufficient  test  time 
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could  enable  the  use  of  standard  rocket  diagnostic  instrumentation  such  as 
load  cells,  thrust  stands,  strain  gauges,  ballistic  pendula,  etc.  Of  course, 
to  obtain  specific  impulse  from  measurements  of  thrust  and  total  impulse  re¬ 
quires  that  the  propellant  mass  flow  rate  and  total  expelled  mass  be  known  ac¬ 
curately. 


Finally,  independent  of  performing  any  new  experiments,  the  existing 
experimental  data  should  be  analyzed  in  greater  detail  through  use  of  the  com¬ 
puter  simulation  code  described  in  Section  2.  Using  the  multiple  pulse  code 
with  time  dependent  laser  absorption  and  real  gas  thermodynamics,  calculations 
should  be  carried  out  for  the  repetitively-pulsed  performance  of  the  conical 
nozzle  operating  with  argon  and  hydrogen  propellants  heated  by  laser  pulse  en¬ 
ergies  and  with  laser  interpulse  times  that  simulate  the  experiments.  It  will 
also  be  extremely  valuable  to  incorporate  into  the  computer  code  a  subroutine 
that  calculates  plasma  radiation  losses. 
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4.  OPTICAL  DESIGN  CONSIDERATIONS  FOR  A  SELF -FOCUSING  NOZZLE 


We  consider  in  this  section  the  optical  performance  of  a  paraboloidal 
self-focusing  nozzle  of  utility  for  laser  propulsion.  The  primary  conclusions 
are: 

•  The  requirement  for  off-axis  performance  allows  the  surface 
to  be  of  only  very  modest  optical  quality. 

•  For  off-axis  angles  £  10°,  only  modest  magnifications  (£25) 
can  be  anticipated. 

•  Optimum  optical  and  fluid  dynamic  performance  appear  to  be 
mutually  exclusive,  consequently  a  compromise  is  required. 

•  The  impact  of  off-axis  focusing  and  breakdown  on  thruster 
performance  and  surface  erosion  should  be  evaluated. 

Details  of  the  analyses  are  given  in  the  following  paragraphs. 

4.1  Aberration  Characterization 

Two  general  classes  of  aberrations  can  be  distinguished:  determinis¬ 
tic  and  random.  The  former  are  those  defects  in  the  image  caused  by  the  de¬ 
parture  of  the  actual  wavefront  from  a  spherical  wave  for  a  mathematically 
perfect  reflective  or  refractive  surface.  The  latter  are  caused  by  departures 
of  the  actual  surface  from  the  desired  mathematical  surface  which  are  random 
because  they  are  associated  with  manufacturing  tolerances.  For  analysis  pur¬ 
poses,  it  is  convenient  to  consider  each  type  separately  because  of  the  dif¬ 
ference  in  their  basic  character. 

4.1.1  Random  Aberrations 

Consider  the  simple  optical  system  shown  in  Fig.  4. 1  which  consists  of 
a  lens  of  aperture  D  and  focal  length  Z.  The  image  of  a  point  object  at  in¬ 
finity  is  given  by 


I(x) 


(Xz)“2  /  dvw(v)U(v)exp[-  “J1  x 


(4.1) 


where  X  is  wavelength,  U(y)  is  the  complex  wave  in  the  aperture  plane  and  W(  jr) 
is  a  simple  aperture  function  (-1  inside  and  0  outside  the  aperture).  We  now 
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assume  that  U(v)  is  characterized  by  the  product  o£  a  deterministic  portion 
(constant  for  a  perfect  aperture)  and  a  random  phase  distortion  («  exp  it(v)). 
The  average  image  plane  intensity  can  be  shown  to  be  equal  to 

2x1  . 

x*A 

<I(x)>  -  p(»)ID(x)  +  — [l-p(«>]  /  dAxD(A)pN(A)e  ~  ~  (4.2) 

(Az) 

where  <...>  represents  a  statistical  average,  ID  is  the  deterministic  image, 
and  E  is  the  energy  collected  by  the  aperture.  td  is  the  optical  transfer 
function  of  the  deterministic  image  (normalized  Fourier  Transform  of  Iq)  and 
PD  is  the  correlation  function  of  the  random  component  of  the  field. 


P(A)  * 

<exp  i[$(v)  -  $(v-A)]> 

p(A)-p(») 

PN(~}  " 

L  1-p(«)  J 

(4.3) 

1  I2  1  I2 

p(«)  - 

|<exp  i4»(v)>  |  =  |*(  1 )  1  , 

where  $  ( 1 )  is  the  first  order  characteristic  function  of  the  random  variable 
♦  (v) . 


In  many  cases  of  interest,  the  correlation  function  is  much  narrower 
than  the  deterministic  transfer  function  leading  to  the  approximate  result. 

)oc 

<I(x)>  «  { I 1 )  1 2  P(x)  +7  [l-|®(  1)  |2]s  (— )}  ,  (4.4) 

(Az)2  *  Z 

when  P(x)  is  the  normalized  point  spread  function  of  the  deterministic  image 
and  S  is  the  spectrum  of  the  random  field  (Fourier  Transform  of  p^). 

From  this  result,  we  note  that  energy  is  conserved  (i.e.,  integral  of 
<I(x)>  over  x  is  equal  to  E) .  Furthermore,  the  energy  in  each  term  is  given 

by 


(Ed/E)  -  I *< 1 > | 2 
(Er/E)  -  1-1 4(1) I2 
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(4.5) 


where  Eq  and  Er  are  the  energies  of  the  deterministic  (first  term  of  Eg.  4.4) 
and  random  (second  term  of  Eg.  4.4}  components  of  the  field,  respectively. 

The  on-axis  intensity  is  given  by: 


<I(o)>  -  -^-r  {l*(1)|2  +  [i-|*(1)|2](A/d)2}  ,  (4.6) 

Uzr 


These  general  characteristics  of  the  image  are  shown  in  Fig.  4.2. 

These  results  indicate  that,  in  principle,  a  deterministic  image  is  always 
present  (of  size  Xz/'D  if  unaberrated).  However,  its  energy  is  attenuated  by 
the  presence  of  the  random  component  of  the  field.  The  random  component  of 
the  field  results  in  an  average  intensity  which  appears  as  though  it  were 
generated  by  a  lens  of  aperture  A,  the  correlation  scale  of  the  random  field. 

The  characteristic  function  depends  on  the  type  of  statistics  involved 
and  the  rms  phase  distortion  (o^).  Values  for  several  common  cases  are  shown 
in  Table  4.1.  As  can  be  seen,  one  wave  (2ir)  of  random  aberration  substantial¬ 
ly  eliminates  the  deterministic  component  in  all  cases.  For  normal  statis¬ 
tics,  a  guarter  and  (1/20)  wave  of  aberration  results  in  a  92%  and  9%  reduc¬ 
tion  in  energy,  respectively.  Conseguently,  high  guality  optical  systems  are 
usually  specified  at  (X/20)  or  better. 

4.1.2  Deterministic  Aberrations 

Except  for  diffraction,  the  on-axis  performance  of  a  mathematical 
paraboloidal  mirror  is  perfect  based  on  third  order  aberration  theory.  For 
off-axis  cases  (i.e.,  direction  of  object  not  parallel  to  axis  of  rotation) 
Coma  and  Astigmatism  occur  (Spherical  Aberration  is  always  zero).  The  angular 
diameters  of  the  blur  circles  associated  with  these  effects  are  given  by: 


Coma:  6^ 


(4.7) 
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TABLE  4.1  First-order  characteristic  functions  for 
various  distributions. 


Distribution 

Characteristic  function 

1  «•(  1 )  1  2 

Value 

2ir 

at  a 

*  4  It 

normal 

2 

exp(-o  ) 

V 

10‘17 

io"69 

uniform* 

2r<,  1/2  , 
sine  [(3)  a,] 

8  x  10~3 

io-4 

4.  .1/2 

-4 

-7 

Simpson 

sine  [(3/2)  '  o.] 

4 

2  x  10 

2  x  10 

2  i 

-2 

-3 

exponential 

<1-0  1 
?> 

2  x  10 

6  x  10 

Laplace 

(1 -0^/2) -2 
$ 

2  x  10  3 

io"4 

gamma 

2  v 

io-2 

00 

1 

o 

Rayleighb 

1F^[1;1/2;-o^/(4-n)] 

10-15 

io’77 

+  ['ff0^/(4-'!r)exp-a^/{4-'rt)] 

y  y 

asincx  =  (sinx/x)* 

blF^(a}  8;  i  X)  «  confluent  hypergeometric  function 


■  I  BJ^nl  where 

n-o 

(cO^Mct+l)  ...  (a+n-1),  (o)q«1. 
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(4.7) 


Astigmatism:  8.  ** 

A 

where  ©  is  the  off-axis  angle,  F  is  the  F-number  of  the  system,  z  is  the  dis¬ 
tance  from  the  paraboloid  apex  to  the  focal  point  and  L  is  the  distance  to  the 
aperture  stop  (see  Fig.  4.3). 


For  the  case  of  interest  here,  a  physical  aperture  stop  is  not  likely 
to  be  used.  Consequently  L  is  defined  by  the  physical  extent  of  the  mirror  or 
the  beam  size.  In  this  case,  the  blur  circle  diameter  due  to  Astigmatism  is 
given  by 


£  r_L_  _  o 

2F  2  ' '  * 

16F 


(4.8) 


The  expressions  given  above  are  reliable  for  modest  apertures  or  fields,  how¬ 
ever,  they  break  down  at  larger  values.  A  very  conservative  estimate  of  the 
resultant  blur  due  to  both  effects  is  given  by  the  sum.  We  will  use  a  root 
sum  square  which  is  perhaps  more  realistic. 


Inspection  of  the  above  forms  lead  to  the  obvious  conclusion  that  for 
off-axis  performance,  a  larger  F-number  is  desirable.  Paraboloids  of  various 
F-number  and  fixed  diameter  are  shown  in  Fig.  4.4.  As  the  F/No.  increases 
with  a  fixed  entrance  diameter,  the  paraboloid  becomes  shallower  and  the  focal 
point  moves  farther  away  from  the  apex.  Above  some  critical  value,  the 
resulting  paraboloid  will  no  longer  function  adequately  as  a  rocket  nozzle. 

The  gas  flow  will  separate  from  the  walls  and/or  the  gas  expansion  rate  will 
not  be  appropriate.  Furthermore,  for  values  in  excess  of  F/0.25,  the  focal 
point  will  lie  outside  the  entrance  aperture.  Unless  some  type  of  mechanical 
extension  is  used,  this  will  likely  yield  an  inefficient  motor  because  the 
thrust  will  not  be  vectored  within  a  small  enough  angle.  Consequently,  nozzle 
design  should  be  based  on  a  trade-off  between  the  conflicting  requirements  for 
optical  vs.  fluid  dynamic  performance. 


The  parabolic  rocket  used  in  laboratory  experiments  (Fig.  4.5)  cor¬ 
responds  to  approximately  F/0.05.  Ray  traces  for  this  thruster  at  0  ■  0°  and 
15°  are  given  in  Figs.  4.6  and  4.7,  respectively.  As  can  be  seen,  the  beam 


71 


5530-2 


4.6  Ray  trace  for  F/0.05  paraboloid  with  0 


does  not  focus  to  any  great  degree  at  6  *  15°.  Ray  traces  for  an  F/0.1  system 
are  shown  in  Fig.  4.8. 


An  approximate  upper  limit  to  off-axis  performance  can  be  obtained  by 
the  ad-hoc  criterion  that  the  aberrated  beam  (rhs  of  Eq.  4.7  and  4.8  converted 
to  spatial  size)  must  be  smaller  than  the  physical  dimension  of  the  parabola, 
i  .e . , 


1M4]2* 

16F 


if?  (±5-’)l2l'/2<4 

16F 


(4.9) 


The  limiting  angle  defined  by  this  expression  is  given  in  Fig.  4.9  as  a  func¬ 
tion  of  F-number.  The  magnification  associated  with  this  spot  size  is  also 
given  in  Fig.  4.9. 


The  results  of  Fig.  4.9  indicate  that  for  low  F-number  systems  ( <  0.1), 
the  laser  beam  will  not  focus  for  off-axis  angles  greater  than  a  few  degrees. 
Optical  performance  at  larger  angles  can  only  be  determined  by  a  detailed  ray 
trace  which  considers  multiple  reflections  or  via  experiment.  Restricting  the 
system  to  small  angles  has  significant  implications  in  terms  of  overall  mission 
configurations. 


4.2  Discussion 

The  analysis  given  in  the  prior  section  indicates  that  because  of 
fluid-dynamic  considerations,  a  low  F/No  paraboloid  will  be  required.  This 
will  result  in  relatively  large  off-axis  deterministic  aberrations  which  will 
produce  a  large  focal  spot  size  at  any  angle  above  a  very  modest  value.  This, 
in  turn,  reduces  the  magnification  dramatically.  Consequently,  any  detailed 
system  study  should  include  the  requirement  that  the  energy  density  of  the 
input  laser  beam  be  within  roughly  a  factor  of  (10-25)  of  that  required  for 
propellant  breakdown. 


When  coupled  with  a  maximum  off-axis  angle  of  10°  or  less,  these  two 
requirements  would  seem  to  be  rather  restrictive  in  terms  of  mission  configura¬ 
tions.  Obviously,  these  requirements  can  be  relaxed  by  implementing  a  larger 
F-number  system.  As  mentioned  previously,  this  has  an  effect  on  performance  as 
a  thruster. 
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Limiting  Angle  (0r) 


5530-8 


Fig.  4.9  Limiting  off  axis  angle  and  corresponding  magnification  for  various 
F/#  paraboloids. 
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Magnification  a  0 . 


To  establish  a  complete  map  of  allowed  operational  parameters  will  re¬ 
quire  a  determination  of  the  limiting  F-number  imposed  by  the  fluid  dynamics. 

While  the  large  spot  size  implied  by  off-axis  performance  limits  the 
achievable  magnification,  there  is  one  advantage  in  terms  of  optical  quality. 
Referring  to  Fig.  4.3,  we  see  that  the  random  image  size  is  given  approximately 
by  (Xz/A)  where  A  is  the  correlation  scale  of  the  surface.  Provided  this 
dimension  is  smaller  than  the  spot  size  associated  with  the  deterministic  off- 
axis  performance,  there  is  no  reason  to  require  a  high  optical  surface  quality 
corresponding  to  a  substantial  fraction  of  energy  in  the  diffraction  image 
(X/10  or  better).  Data  from  one  vendor  of  paraboloidal  mirrors  suggests  that 
this  scale  is  of  order  0.1  mm  which  is  consistent  with  other  data  relative  to 
surface  ( non-optical )  finishes.  For  a  one  meter,  F/0.05  thruster,  this  value 
yields  a  random  image  size  of  less  than  1  mm  at  a  wavelength  of  0.35  y  which  is 
much  smaller  than  the  20  cm  nozzle  size  at  nominal  focus  (for  a  one  meter  aper¬ 
ture).  Consequently,  several  waves  of  random  aberration  could  be  tolerated 
without  serious  impact  on  performance.  To  determine  the  actual  amount  allowed 
would  require  an  alternate  analysis  which  takes  into  account  surface  slope 
(i.e.,  0^,/A ) ,  and  shadowing  effects  which  can  become  important  for  low  inci¬ 
dence  angles  on  a  surface  with  steep  slopes.  The  analysis  of  Section  4.1  does 
not  adequately  account  for  these  effects,  i.e.,  it  is  appropriate  for  a  moder¬ 
ately  good  surface.  This  alternate  analysis  is  available  in  several  standard 
texts  but  it  is  probably  not  necessary  to  pursue  this  issue  in  detail. 

One  effect  which  has  not  been  evaluated  is  the  impact  of  off-axis 
focusing  and  breakdown  on  thruster  performance.  As  the  beam  moves  off-axis, 
the  focal  point  also  moves  off-axis.  This  may  be  important.  However,  as  in¬ 
dicated  in  the  prior  results,  the  beam  quickly  becomes  large  and  so  it  may  not 
be  appropriate  to  conceptualize  the  situation  in  terms  of  breakdown  occurring 
at  a  small  point  which  is  located  at  a  substantial  distance  from  the  axis. 

A  related  issue  is  surface  erosion  which  may  occur  as  the  focus  and 
breakdown  region  moves  close  to  the  walls.  Again,  however,  it  is  probably  in¬ 
appropriate  to  consider  the  problem  in  this  fashion.  Rather  a  more  conserva¬ 
tive  approach  would  be  to  determine  the  effect  on  the  surface  when  the  break¬ 
down  occurs  in  a  large  volume  which  has  a  dimension  of  order  of  the  nozzle  in 
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the  vicinity  o£  nominal  focus.  It  seems  reasonable  that  this  assumption  would 
lead  to  a  reasonable  estimate  of  surface  erosion  effects. 

In  summary,  it  appears  that  the  allowed  off-axis  angle  is  significant¬ 
ly  limited  by  the  low  F-numbers  which  probably  are  required  for  suitable  fluid 
dynamic  performance.  This  has  the  further  effect  of  yielding  large  spot  sizes 
which  limit  the  achievable  magnification.  This  does  have  benefit  in  that  only 
very  modest  optical  surface  quality  is  required.  To  establish  a  complete  map 
of  allowed  operational  conditions  will  require  the  evaluation  of  the  maximum 
allowed  F-number  which  is  controlled  by  the  fluid  mechanics.  Off-axis  focus¬ 
ing  could  be  a  problem,  and  a  conservative  approach  to  evaluating  surface 
erosion  would  be  to  assume  that  breakdown  occurs  in  a  volume  with  dimension  of 
order  of  the  nozzle  size  in  the  vicinity  of  nominal  focus.  It  should  be  noted 
that  these  conclusions  are  based  on  a  rather  simple  analysis.  Depending  on  the 
input  beam  profile,  adequate  energy  densities  which  will  lead  to  gas  breakdown 
may  exist  in  a  small  region  even  for  large  angles  and  nominal  spot  sizes.  This 
effect  could  be  significant  and  hence,  a  detailed  ray  trace  which  includes  beam 
profiles  should  be  implemented  to  provide  a  complete  assessment  of  these 
issues • 


